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Recent data from angle-resolved photoemission experiments published by Zhou et al. [Nature, 
423, 398 (2003)] concerning a number of hole-doped copper-oxide based high-temperature supercon- 
ductors reveal that in the nodal directions of the underlying square Brillouin zones (i.e. the directions 
along which the d-wave superconducting gap is vanishing) the Fermi velocities for some finite range 
of k inside the Fermi sea and away from the nodal Fermi wave vector fep are to within an experimen- 
tal uncertainty of approximately 20% the same both in all the compounds investigated and over a 
wide range of doping concentrations and that, in line with earlier experimental observations, at some 
characteristic wave vector fc* away from fep (||fe* — k^W typically amounts to approximately 5% of 
llfcrll) the Fermi velocities undergo a sudden change, with this change (roughly speaking, an increase 
for ||fe|| < ||fe*||) being the greatest (smallest) in the case of underdoped (overdoped) compounds. 
We demonstrate that these observations establish four essential facts: firstly, the ground-state mo- 
mentum distribution function n(fe) must be discontinuous at fe = fe*; with and v^^ denoting 

the measured velocities close to fc = fc*, |jfc^ || < ||fe*|| < ||fej||, and = n{k^) — n(fc+) > 0, we 
obtain 

n(K^. ) 

in which is the 'Fermi' velocity corresponding to the case in which = (for two- 

body interaction potentials of shorter range than the Coulomb potential, v/^^ ~ v^^ whereby 

< v'^_ ^21)^^); secondly, the single-particle spectral function A{k;e) must at fc = fe* pos- 
sess a coherent contribution corresponding to a well-defined quasiparticle excitation at an energy 
of approximately 70 meV below the Fermi energy; thirdly, the amount of discontinuity Zj^^ in n(fc) 
at the nodal Fermi points must be small and ideally vanishing; fourthly, the long range of the 
two-body Coulomb potential is of vital importance for realization of a certain aspect of the ob- 
served behaviour in the Fermi velocities (specifically) in the underdoped regime. The condition 
^fcp = conforms with the observation through an earlier angle-resolved photoemission experiment 
by Valla et al. [Science, 285, 2110 (1999)], on the optimally doped compound Bi2Sr2CaCu208+i, 
which shows that the imaginary part of the self-energy along the nodal directions of the Brillouin 
zone and in the vicinity of the nodal Fermi points satisfies the scaling behaviour characteristic of 
marginal Fermi-liquid metallic states, for which Z/^^ is indeed vanishing. We present arguments 
advocating the viewpoint that the observed 'kink' in the measured energy dispersions cannot be a 
direct consequence of electron-phonon interaction, although a finite Zf,^ may possibly arise from 
this interaction. In other words, even though possibly vital, the role played by phonons in bringing 
about the latter 'kink' must be indirect. Our approach further provides a consistent interpretation 
of the observed sudden decrease in the width of the so-called 'momentum distribution curve' on ||fe|| 
increasing above ||fe*||. 
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§ 1. Introduction 

Angle-resolved photoemission spectroscopy (ARPES) 
[1,2] provides information with regard to the single- 
particle spectral function Ac,{k;e) of systems for £ < /x, 
where e is the external energy parameter, ^ the chem- 
ical potential, fe the wave vector and a the spin index. 
The quantity measured experimentally is the photoelec- 
tron intensity Ia{k; e) which at finite temperature T and 
under some specified approximations is proportional to 
/(£-;u)A<,(fe;£), where /(e) = l/(l-Fexp(£/[A;Br])) is the 
Fermi distribution function [1,2]. The experimental re- 
sults obtained by Zhou et al. [3] and other workers [4-8], 
to be discussed in the present paper, are deduced from 
the so-called 'momentum distribution curve' (MDC); this 
curve is obtained by setting out Ia{k;e) along the de- 
sired directions in the fc space for a fixed value of e. 
A related curve, referred to as the 'energy distribution 
curve' (EDC), is obtained by plotting Ia{k; s) for a fixed 
fc along the e axis. Amongst others, the determination 
of the peak positions in Acr{k;e) from the correspond- 
ing MDC does not suffer from the ambiguities associated 
with /(c - /i). 

In this paper we first consider the relationship between 
the energy dispersions as measured in angle-resolved pho- 
toemission experiments and the single-particle excitation 



energy dispersions as determined by solving the quasi- 
particle equation in terms of the self-energy operator; the 
latter accounts for the correlation in the ground states 
(GSs) of the underlying systems. The formalism concern- 
ing the measured energy dispersions that follows from our 
purely heuristic arguments turns out to coincide, in all 
details, with that developed by the present author in his 
recent investigations regarding the nature of the uniform 
metallic GSs of the conventional single-band Hubbard 
Hamiltonian [9] and more general single-band Hamilto- 
nians [10] in which the particle-particle interaction can 
be of arbitrary range. The specific aspects of the formal- 
ism developed in [9,10] that greatly suited the aims in the 
aforementioned investigations are that for the indicated 
GSs it formally yields 

(i) the exact GS total energy, 

(ii) the exact Fermi energy, 

(iii) the exact Fermi surface and 

(iv) the exact momentum distribution function. 

The formalism has made possible [9,10] to make exact 
predictions concerning the values that the GS momen- 
tum distribution function na-{k) can take for k in the close 
neighbourhoods of the underlying Fermi surface >Sf;<7- On 
the basis of this, it has been possible to expose [10] the 
significant influence that the range of the two-body inter- 
action potential can have on na-{k) and consequently on 
the partitioning of the GS total energy into 'kinetic' and 
'interaction' parts. This observation led us to the conclu- 
sion [10] that the experimentally-inferred excess 'kinetic' 
energy in the normal states of the copper-oxide based 
high-temperature superconductors [11,12] (for other per- 
tinent references see [10]) may be signifying the impor- 
tance of the long range of the two-body interaction poten- 
tial in determining other distinctive aspects of the nor- 
mal metallic states of these compounds that mark them 
as unconventional metals [13]. 

For the normal uniform metallic GSs of the single-band 
Hubbard Hamiltonian, we have established that, so long 
as these are Fermi liquids, for fe inside the Fermi sea and 
close to the Fermi surface, a single-particle energy dis- 
persion £^.^, which in [9,10] was attributed to some 'fic- 
titious' particles and in this paper we reason to be the one 
measured through the ARPES, is up to a multiplicative 
fc-independent constant equal to the energy dispersion 
efe;o- of the Landau quasi-particles inside the Fermi sea 
(more precisely, the energy dispersion e^.^ as obtained 
from the quasi-particle equation in Eq. (3) below; see 
§ 3.1); similarly for e^.^ (Eq. (23)) and e'^.^ concerning fc 
outside the Fermi sea. For cases where the two-particle 
interaction potential is the long-range Coulomb poten- 
tial, in [10] it was however shown that no such strict rela- 
tionship exists between e^.^ (and similarly £^.^) and e^.^ 
(^fe-o-)' eyen for Fermi liquid metallic states. In this 
paper we rigorously demonstrate that in spite of these as- 
pects, in cases where for a given fe = fe* the quasi-particle 
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equation (Eq. (3) below) has solutions and 
corresponding to well-defined Landau quasi-particles, ir- 
respective of the nature of the two-body interaction po- 
tential identically coincides with e'^ (if e^.^ is 
continuous at A; = fe*; Eq. (27)) and with 
(if e^.^ is continuous at fc = fe*; Eq. (30)); we present 
the appropriate expressions for and ^Xi.-cj terms 
of e^j, and e^,: respcctivclv for cases in which either 

or both of e^.^ and e^.^ arc discontinuous at fc = A;*; 
according to our analysis (§ 3.4), these functions cannot 
both be continuous at k = ki,. 

In the light of the above-mentioned similarities be- 
tween £^.^ and the possible solutions e^.^ of the quasi- 
particle equation associated with particle-Uke excitations 
(as opposed to collective low-lying excitations in non- 
Fermi-liquid metallic states) as well as the heuristic argu- 
ments underlying the definition of e^.^ (see § 2), in this 
paper we assert that the ARPES data concerning the 
dispersion of the single-particle excitation energies are 
to be compared with e^.^. The fact that the behaviour 
of e^.^ corresponding to Fermi-liquid metallic states of 
systems of fcrmions interacting through the long-range 
Coulomb potential in the neighbourhood of S-^-^a demon- 
strably differs from that of the asymptotic solution (sec 
§ 3.1) e^.^ of the quasi-particle equation in the neighbour- 
hood of <Sf;o- (see above) does not diminish the strength 
of our arguments. In this connection we should empha- 
size that our assertion has bearing on the energy dis- 
persions as measured hy means of the ARPES and thus 
does not negate the fundamental significance of the pos- 
sible asymptotic solutions of the quasi-particle equation, 
nor the essential role that these (explicitly, the asymp- 
totic solutions corresponding to the exact quasi-particle 
solutions for fc G Sp-a] § 3.1) play in determining the low- 
temperature thermodynamic and transport properties of 
the Fermi-liquid metallic systems [14]. 

Our point of departure in this paper consists of the con- 
sideration that the measured (by means of the ARPES) 
single-particle energy dispersions coincide with the ex- 
pectation value of e with respect to the normalized distri- 
bution of the single-particle excitations whose energies lie 
above —E and below fi, that is A„{k; e)/ J^^ de A„{k; e) 
(note the dependence on fc of the normalization ampli- 
tude); for systems with bounded single-particle energy 
spectrum (such as those described by the single-band 
Hubbard and the extended Hubbard Hamiltonians), E 
can be identified with infinity (§ 2). This considera- 
tion regarding the measured energy dispersions provides 
the opportunity to introduce a quantitative size for the 
width of the main peak in the single-particle spectral 
function Arr{k:e) centred at e = e^.^, which is express- 
ible in terms of Ji^^de Acr{k;e), J^^de s"^ Acr{k;e) and 
e^.g. (see Eq. (7)). According to this measure and under 
the conditions that we shall explicitly specify in the due 



place, the behaviour of the width of the single-particle 
spectral function centred at e^.^ is directly correlated 
with the behaviour of n^{k); on transposing fc through a 
region where na-{k) undergoes a rapid decrease (increase), 
this behaviour is directly reflected in a concomitant de- 
crease (increase) in the latter width. This, which for 
metallic states can be shown to be the characteristic as- 
pect of the corresponding A^-^k; s) for fc approaching the 
underlying 5f;<7, ^ thus turns out to be a generic property 
of Aa-{k:e) pertaining to systems of interacting fcrmions 
for fc in the neighbourhoods of all points (and not solely 
the points of S^-a) where n^{k) exhibits large local vari- 
ation. According to this insight, the experimental ob- 
servations in [3,5,7,8] with regard to a sharp decrease in 
the width of the peak in A^{k; e) centred at e = £^.^ for 
fc = fci, is indicative that fc^, is a point at which the un- 
derlying GS momentum distribution functions undergo 
a sharp decrease, if not discontinuity. Our analysis of 
other aspects of the experimental observations reported 
in [3-8] supports a discontinuity in the underlying no-(fc), 
with n„{k^) > no-(fe^); where ||fe^ — fe^|| is infinites- 
imally small and fe^ denotes the vector nearest to the 
nodal Fermi wave vector kp-a-- 

In the sections that follow, we combine our theoreti- 
cal considerations with brief discussions of the pertinent 
experimental observations in [3-8] . Our theoretical treat- 
ment in this paper is restricted to T = 0. 

§ 2. Preliminaries 

Assuming that, for a given fe, the experimentally deter- 
mined Aa{k;e) is dominantly sharp at and symmetrical 
around e = Sk-.a (later we shall identify this energy with 
e^.^ for reasons that will be clarified), Sk-.a would coin- 
cide with the expectation value of e with respect to the 
energy distribution function t]r^Acr{k;e). In this con- 
nection it is important to realize that A^{k;e) is posi- 
tive semidefinitc, as befits a proper distribution function; 
however, although h^^ Aa-{k; e) is normalized to imity 
for s over (— oo, oo), this is not the case for s restricted 
to the semi-infinite interval {—oo, fj]. Thus, under the 
above-mentioned conditions, for Sk-a < corresponding 
to 'occupied' single-particle states, one should have (see 
Eq. (39) below) 



^ On the basis of experimental observations, it has even been 
advocated [15,16] (see however the cautionary remarks in ref- 
erence 8 of [16]) that Fermi surface be defined as the locus of 
points in the k space at which )|VfenCT(fe))| is ma^ximal. This 
definition has been utilized in [17] where the n^ik) pertaining 
to the GS of the t-J Hamiltonian has been calculated from a 
high- temperature series expansion. See also [9]. 
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where 



1 



(1) 



(2) 



stands for the GS momentum distribution function per- 
taining to particles with spin index a. The right-hand 
side (RHS) of Eq. (1) identically coincides with the ex- 
pression for the energy dispersion e^.^ as defined in ^ 
[9,10]. We note in passing that e^.^ < for all k [9,10]. 
Evidently, the experimentally measured Aa{k: s) exhibits 
more structure than a single dominant peak, however so 
long as these additional structure correspond to high- 
energy processes (e.g. core-level peaks, which is not de- 
scribed by single-band models), for a relatively small 
range of k, such structure can only account for a rigid 
(i.e. fe-independent) shift of the RHS of Eq. (1) with 
respect to the low-energy peak position of the experi- 
mentally determined Aa{k;e). 

Although e^.^ (i.e., the RHS of Eq. (1)) is not formally 
a 'quasiparticle' energy, in the sense of unconditionally 
being the solution of the quasiparticle equation (see later) 



Ek + hT,^{k;e) = e. 



(3) 



in which Sk is the single-particle energy dispersion de- 
scribing the non-interacting system (see Eq. (36) below) 
and E(j(A;;e) the self-energy, in [9,10] it has been rigor- 
ously shown ^ that, for k approaching the interacting 
Fermi surface S^-cr, s^-a approaches the exact Fermi en- 
ergy Ef; with fcpicr e S^-a and ^ fc^^ {kt-a) inside (outside) 
the Fermi sea FS^ (in this paper we denote the set of k 
points complementary to FScr with respect to the avail- 
able wave-vector space by FSo- so that k^.^ e FScr) and 
infinitesimally close to k^-a, we have [9,10] 



Ef = 



j:!^deeA,{kf.y,e) 
j!l^deA^{kf^^;e)- 



(4) 



The freedom in the choice of either of the two vectors k^.^ 
and fcjto- in the above expression reflects the continuity 
of e^.^ in a neighbourhood of Sp-a (for details see §§ 3.4 
and 4.2 below); our use of kf.^, rather than simply k^-a 
(the continuity of s^.^ at fc = k^-^ notwithstanding), 



^ See, for instance, Eqs. (6) and (39) in [10]. 
^ See, e.g., Eq. (40) in [10]. 

* For reasons that we have presented earlier [18,9], one can 
identify k^.^ with kF;a (i.e., SF;a can be considered as a proper 
subset of FSct); on the other hand, strict distinction has to be 
made between kt^a and fcF;^- The same applies to fe* that we 
frequently encounter in the text. See §§ 3.1 and 3.2. 



is prompted by the fact that, in cases where ncr(fc) is 
discontinuous at fe = k^-a, the expression on the RHS of 
Eq. (4) is in its present form formally ill defined. 

We point out that for metallic GSs, e = Ef satisfies 
Eq. (3); Fermi surface Sp-a is the locus of all k points for 
which this equation is fulfilled for e = £f and moreover 
Sk + ?iScr(A;;eF) — has different signs for k infinitesi- 
mally transposed in opposite directions along the normal 
to the constant-energy surface, i.e. S^-cr, described by 
Sk + ?iScr(fe;£F) = £f- ^ Equation (4) presages the pos- 
sibility that the RHS of Eq. (1) may at least be a useful 
mathematical device for the purpose of estimating the 
energies of the low- lying single-particle excitations in in- 
teracting metallic systems. 

As a measure both of the accuracy with which e^.^ de- 
scribes the location of the main peak in A^rik; e), e < n, 
and of the width of the latter peak, one can consider the 
variance Ae^.^ of s. Since for the systems in which Sk 
is unbounded from above, which for the uniform GSs 
considered in this paper is feasible only for those de- 
fined on the continuum (to be contrasted with those de- 
fined on lattice), one has Aa{k;e) ^ ?i^So-(fc)/(7r^|e]^) 
for |e| oo (see Eqs. (239a) and (239b) in [19]), where 
S(^(fe) stands for a well-defined function, ^ it follows that 
for these systems /^^ de A„{k; e) is logarithmically di- 
vergent for E oo. Consequently, for systems in which 
Efc is unbounded from above it is necessary to employ a fi- 
nite energy cutoff, E, in order for Ae^.^ to be a meaning- 
ful quantity. ^ Thus, although for the latter systems the 
average value of e, that is £^.^, is bounded for £^ — > oo, 
a consistent formalism for the determination of the posi- 
tion of the fundamental peak and the associated width in 
the single-particle spectral function requires introduction 
of a single cut-off energy throughout. We thus define (see 
Eqs. (1) and (39)) 



h ncr{k; E) 
where (see Eq. (2)) 

1 

n„{k; E):=- J de Aa{k;e). 



(5) 



(6) 



® For completeness, for a given Sf;ct there also exists a set 
(S^o^, a generic point of which we denote by fc^j, (see footnote 
4), for which Eq. (3) is satisfied at e = e^, with infinitesi- 
mally greater than Ef. See §§ 3.1 and 3.2. 

® See Eq. (227) in [19] and Eq. (A5) in [10]. By the as- 
sumption of the stability of the underlying GSs, Sa{k) is 
non-negative. 

^ It is evident that, by defining Ae^._^ as the expectation 
value of, for instance, (e — e^.^l, no need for the introduction 
of an energy cut-off would arise, this is, however, at the ex- 
pense of forsaking the simplicity of the expression for Ae^.^ 
corresponding to the standard variance of e. 
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For the variance we consequently have 

v2 



1/2 



1/2 



For systems with bounded Sk, E in the above expressions 
can be equated with +00. For systems with unbounded 
Efc, it is necessary to set E equal to some finite value; 
however, none of the main results in this paper crucially 
depends on the precise value chosen for E so long as this 
value is greater than the absolute values of the energies of 
interest, specifically t^.^. Consequently, in the remaining 
parts of this paper we leave out E as the argument of the 
pertinent functions; for consistency with the results in 
[9,10] (which do not consider Ae^.^), at places wc still 
write —00, rather than —E, in the lower boundaries of 
the integrals over e. 

It is interesting to note that according to the RHS 
of Eq. (7) a smooth (or continuous) behaviour in 
J^^de e'^A„{k;e) and e^.^ in the neighbourhood of a 
point fey,, where n^{k) undergoes a steep (or discontinu- 
ous) change, must be accompanied by a decrease or in- 
crease in Ae^.^ (corresponding to a narrowing or broad- 
ening of the peak centred around £k^,-a) upon transpos- 
ing k from k~ to k^ , depending on whether na{k~) > 
"cr{kt) or r\cr{k~) < n^{k^) respectively. This observa- 
tion is based on the non-negativity of the argument of 
the square root function on the RHS of Eq. (7) for all 
k. In fact, in § 3.7 we rigorously demonstrate that, at 
points where e^.^ is continuous (such as the points of 
the underlying Fermi surface; see § 3.4) and ncr(fc) is dis- 
continuous, Ae^.^ must necessarily be discontinuous. In 
the light of this, the sudden decrease in the width of the 
MDC at fe = fe* (while transposing k towards the nodal 
Fermi point kp.^r along the nodal direction of the first 
Brillouin zone (IBZ)) in the experimental observations 
by Zhou ct al. [3] (corresponding to the hole concentra- 
tion X — 0.063) is indicative of a sharp decrease (if not 
discontinuity) in r\a{k) at fe = fc*. The same applies for 
the similar observations in [5,7,8]. 

It is important to realize that, for interacting GSs, 
docs not vanish for any k, not even for k — > Sf-^o-- 
This follows from the fact that, with the exception of 
e = £^.^, (e — e^.^)^ is positive everywhere along the 
e axis and that Aa{k;e) > 0, Ve. For non-interacting 
GSs, where Ac,{k; e) consists solely of a single 5 function 
(see appendix A), the corresponding Ae^.^ is naturally 
identically vanishing. 

In order to appreciate the significance of Ae^.^, it is ap- 
propriate to consider the following. Fermi-liquid metallic 
states are characterized by smooth quasiparticlc energy 
dispersions Sk-.a in some neighbourhood of the underly- 
ing Sp-cr (see § 3.1), a property that on general grounds 



one expects to be shared by e^.^ (see § 3.4). For con- 
ventional Fermi-liquids, S^-^a- constitutes the locus of the 
k points in the vicinity of which the underlying n„{k) 
undergoes its strongest variation (see footnote 1). With 
these in mind, our above-mentioned statement with re- 
gard to the behaviour of Ae^.^ is seen to be in full con- 
formity with the fact that for the Fermi-liquid metal- 
lic states the width of the quasi-particle peak, centred 
at e = £fe;CT, in the corresponding Acr{k;e) also dimin- 
ishes. It is important at this point to realize that the 
descend towards zero of |Im[SCT(fc; e)] | (which in the case 
of the conventional Fermi liquids is like (sk-cr — £f)^) for 
k approaching Sp-^a- does not imply the same for Ae^.^ 
(which otherwise would contradict our above statement 
with regard to the strict positivity of the latter quantity 
in the case of interacting GSs), for Ae^.^ takes account 
of Ac{k;£) in its entirety and not solely of the coherent 
contribution to A„{k;e) (i.e. S~{s) in Eq. (8) below) 
which indeed has a vanishing width. 

In the light of the above observations, we conclude 
that the standard variance Ae^.^ defined in Eq. (7) pro- 
vides a reliable representation of the width of the peak 
in the single-particle spectral function Acr{k;£) centred 
at £ = £^.^ for k close to any point (which may or may 
not be a point of Sp-^a) at which e^.^ is continuous and 
the underlying n„{k) undergoes a sharp or discontinuous 
change. 



§ 3. Generalities 
3.1. Some remarks on the single-particle excitation energies 



For an arbitrary wave vector k the quasi-particle equa- 
tion, Eq. (3), has in general no solution. ^ In spite 
of this, in cases where k is located in a small neigh- 
bourhood of a point, say fco, where Eq. (3) is satisfied 
and the self-energy is sufficiently smooth with respect 
to k (for e = eka;a) and £ (for k = ko), an energy 
dispersion ek;a is capable of being defined through an 
asymptotic series in terms of the asymptotic sequence 
{{k — feo)^ IK — 0,1,...}, with the zeroth-order term 
in this series coinciding with the exact solution Sko;a of 
Eq. (3); we thus refer to £k;a as the asymptotic solution 



* We have in various papers elaborated on this subject mat- 
ter, for which we refer the reader to [19] (§ 3.4) and the ref- 
erences therein. Here and in the following when stating that, 
for an arbitrary k, in general, Eq. (3) does not possess a solu- 
tion, we do not thereby consider the possible solutions of this 
equation on the non-physical Riemann sheets whose determi- 
nation requires the knowledge of the analytically continued 
Y:a{k;e) into these Riemann sheets (see, e.g., [20]). 
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of Eq. (3), reflecting the fact that, for the k at issue, 
Eq. (3) has no true solution. By the non-analytic nature 
of J^cr{k;£) at (fe,e) = (feo^feoit^) (^^^ later), the afore- 
mentioned series can only be of flnitc order, that is there 
exists some finite integer £q such that the coefficient of 
the ^th-order term in this series will be unbounded for 
^ > ^0- Fermi-liquid metallic states arc those for which 
the £o corresponding to all k in the neighbourhood of 
the underlying Fermi surface S^-cr is strictly greater than 
unity [18,19,9]; the effective mass, or the Fermi velocity, 
of the Landau 'quasi-particle' at k is related to the co- 
efficient of the asymptotic term linear in (fe — feo), with 
feo S Sp-a', for anisotropic metallic GSs the latter coeffi- 
cient is in general a non-diagonal Cartesian tensor. 

Concerning the above-mentioned non-analytic nature 
of Eo-(fe;e) at (fc,e) = (ko, Sko-.a), this follows from the 
fact that for interacting systems (specifically those in the 
thermodynamic limit) the solutions of the quasi-particle 
equation are not isolated points but accumulation points 
[18,19]; this aspect is directly related to the overcom- 
pleteness of the set of amplitudes of the single-particle 
excitations (the so-called Lehmann amplitudes) in the 
single-particle Hilbert space of interacting Hamiltonians 
(see [19]). 

In view of the specific applications in this paper, we 
mention that even for single-band models, based on a 
non-interacting energy dispersion (in this paper denoted 
by Efc), the possible solutions of Eq. (3) for a given fc oc- 
cur in pairs (for a more precise specification sec footnote 
9 below), of which one is associated with point k~ and 
one with point fe"*" where k~ and k^ are infinitesimally 
close to fe (see Fig. 1). This aspect, which is universally 
neglected in the literature for k G Sp-^ (see however [18] 
and specifically [9]; sec also footnote 5 above), attains 
an extraordinary significance in cases where the above- 
mentioned fe is located at some finite distance from S^-a 
(considering for the moment metallic GSs), and hence 
our explicit reference to ej.^ = s^T-a in § 1 (see Eq. (13) 
below). This aspect will be fully clarified in § 3.2 below. 

3.2. Regular and singular contributions to the single-particle 
spectral function Aa{k;s) 

In anticipation of our following discussions and in view 

of Eq. (4) it is instructive to introduce a decomposition of 
Aa-{k;£) in terms of its regular ('incoherent') and singu- 
lar ('coherent') parts. To this end, let fe* denote a wave 
vector at which na-{k) is singidar (not necessarily dis- 
continuous). For fc = fc^, where fe^ are infinitesimally 
close to ki, and whose distinction is implicitly determined 
through Eq. (9) below and the accompanying conditions 
£f^T.^ > A*' '^''^^ write 

A,{kf;e) = R„{e) + S^{e), (8) 

where Rcr{£) {S^{s)) stands for a regular (singular) 'func- 
tion' of s; for fe* coinciding with a regular point of ncr(fc). 



SJ(e) would be identically vanishing. We point out that 
the equality in Eq. (8) is between two distributions (as 
opposed to functions) so that it applies in some space of 
test functions. In the event that ncr(fc) is discontinuous 
at fc*, in the simplest case ^ one has 

S^ie) = hZ^^.J{e-e^^J, (9) 

where Zj^^f.^ > 0, e^T.^ > M i^oi the latter see the text 
following Eq. (11) below), and 

The strict equality on the left-hand side (LHS) of Eq. (10) 
(whereby in the subsequent parts of this paper we denote 
■^fe^F-CT follows from the normalization condition 

1 f°° 

-J dsA^{k;e) = l, (11) 

which applies for all fc, including fc~ and fe^. We point 
out that £^--0- ^^'^ ^k+-a satisfy Eq. (3) for k — k^ and 
k = k~l respectively (see Fig. 1) and, moreover, the ex- 
pression for S;]r(e) in Eq. (9) is specific to cases where 
T,„{kf; £) is continuously differentiable with respect to £ 
in the neighbourhood of £ = e^t.,^ (c/. Eq. (80) below); 
see [18]. 

The normalization condition in Eq. (11) in combina- 
tion with the defining equation for rtcr{k) in Eq. (2) imply 
that 

1 f°° 

-J d£A<,(fe;£) = l-n^(fe). (12) 

As is evident, the discontinuity of ncr(fe) at fe = fc* re- 
quires that £k--a- < A*' for otherwise S^(e) would not 
have any impact on the behaviour of nCT(fe), defined ac- 
cording to Eq. (2). On the other hand, the discontinuity 
of na{k) at fe = fe* implies discontinuity of 1 — r\(^{k) at 
k = ki,, which, following Eq. (12), can only materialize 
if the condition £k+-a- > M is satisfied [9]. 

Since for fe* located at some finite (i.e. non-vanishing) 
distance from Sp-^a the solutions £j,-.g. and £k+-a differ 
by some finite amount ^° (see §§ 4.2 and 4.3 as well as 



^ In general, one may have SJ(£) = ^iXljJi ^(^ " 

e^^^ ) where, for the same reasons that -^j.-.^ = -^fe+.CT 
(Eq. (10)), {Z^^l } and {Z^^l } are interdependent. As con- 

sidorations involving AfJ^ > 1 only unnecessarily complicate 
the analysis, in this paper we explicitly deal with cases where 
< 1. 

Since by the arguments presented in the text following 
Eq. (12) we have £k~ a < M < ^fc+ a' follows that in the 
event that e^^- and e^^- would differ only by an infinitesi- 
mal amount, fe* would either be a point of .Sp;^ or infinitesi- 
mally close to it. 
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Fig. 1), it is of crucial importance that particular atten- 
tion is given to the superscripts — and + in s^t.„ which as 
a matter of habit one is apt to ignore. We shall elaborate 

on the single-particle energies s^T-cr and the relationship 

< 

between these and e^.^ in § 3.5 where we are equipped 
with the details necessary for this task. Suffice it for the 
moment to mention that the single-particle excitation en- 
ergies of single-band models can be viewed as consisting 
of at least (see footnote 9) two branches, conveniently 
denoted by Si^t-„ with Sk-.^ < fJ, and Si^+.„ > /U. The 

energy dispersions e^.^ dealt with in this paper are 'cari- 
catures' of e^T-a which, however, when continuous at the 
points of discontinuity of n„{k) (see § 3.4), coincide with 
^kT-a (c/- Eqs. (27) and (30) below). From this perspec- 
tive, it is more transparent to denote Sj^^:.,j by s^.^- In 
this paper we shall therefore alternately use the following 
convention: 



'fc;(T — 



(13) 



We emphasize that, for a general k, Eq. (3) has no solu- 
tion so that, in general efeT;^, or ej.^, should be consid- 
ered in the manner specified in § 3.1. 

From Eqs. (8) and (9) and the condition e^^.^^ IJ- 
which implies S^(£) = for £ > /i, for cases where no-(fe) 
is discontinuous at fe = fe* one obtains 



R„{e) = A„{k^\e) for e</i. 



(14) 



Prom this result as specialized to fe* = kj,.a e S^-cr and 
Eq. (4), for interacting systems we immediately obtain 



£f = 



/^^d££R.(£) 



for ki, G S^-a- 



(15) 



This result is of experimental significance, as the pos- 
sible singular contribution to A„{k'f.^;e) is not capable 
of being resolved in real experiments. It is also of sig- 
nificance to our subsequent considerations as it shows 
that, for k sufficiently close to S^-,^, there is no restric- 
tion on the actual shape of the corresponding R„{e) in 
order for the expression on the RHS of Eq. (1) to amount 
to a reliable description of the dispersion of the low-lying 
single-particle excitations in interacting systems, this in 
surprising contrast with the heuristic arguments in § 2 
that initially motivated our introduction of the expres- 
sion in Eq. (1). For illustration, assuming that no-(fc) is 
discontinuous at fe* = k^-a € <5f;<7 and denoting Zk^.„ by 
Zk^_^, from Eqs. (2), (8) and (9) we obtain the following 
relationship: 



With reference to our remark in footnote 9, we note that 
for cases where S~ (e) has the more general form as presented 



Zk,.„ = n.(fe-J - n<,(fe+J, (16) 
which is a well-known Migdal [21,22] theorem. 

3.3. On ratios of discontinuous functions 

A result that proves considerably useful in our consid- 
erations in this paper is based on some simple observa- 
tions. Let 



h{x):= 



(17) 



and assume that g{x) is discontinuous at a; = a;o so that 
Ag ^ 0, where 



(18) 



in which Xq and Xq are infinitesimally close to xq G 
{xq ,Xq). Assuming h{x) to be continuous at a; = xq, i.e. 
h{x^) = h{xo), one readily deduces that 



where 



h{xo), 



A/:=/(a;o ) - /(xo+). 



(19) 



(20) 



Note that, according to Eq. (19), milcss h{xo) — 0, dis- 
continuity of g{x^ at a; = a:;o implies that of /(a;) at 
X = .To, that is Ay ^ 0. A mmibcr of applications of the 
result in Eq. (19), with remarkable consequences, will be 
encountered in the following sections. 

Another result that will prove useful to our later anal- 
ysis concerns a situation in which h{x) as defined in 
Eq. (17) is discontinuous at a; = a;o, similar to f{x) and 
^(a;). Let 

^h:=Kx-^)-h{xl). (21) 
Through some straightforward algebra one obtains 

A/-Ma^o)Ag 



Aft 



(22) 



One observes that the requirement A/, = immediately 
leads to the result in Eq. (19). 



in that footnote, the counterpart of the result in Eq. (16) 
would in principle differ from that which stands in Eq. (16). 
However, the considerations in § 3.6 demonstrate that the 
corresponding to fc* G iSf;o- can at most be equal to unity. In 
other words, Eq. (16) is more general than the simple form for 
S~(e) in Eq. (9) would suggest. For completeness, in Eq. (16), 
]!if.y.=k-e;rj^ i^nilK.-e;!,) 1 K i 0, where nift-e^a) is the unit vector 
normal to (Sp;^ at fe = feF;^, pointing from the inside to outside 
the Fermi sea. 



7 



3.4. Possible discontinuities in ef: and 

Some of the crucial results in this paper rely on the 
condition of continuity of e^.^ at such points as fe = 
fe* where ncr(fe) is discontinuous. Here we investigate 
whether the two conditions can be compatible. To this 
end we introduce e^.^, defined as (see Eq. (48) below) 

It can be easily shown that e^.^ and s^.^ satisfy the fol- 
lowing exact relationship [9,10] 

"a{k)e<^ + {l-n„{k))s>^=ek + hK'{k), Vfc, (24) 

where E^''(fc) stands for the Hartrcc-Fock self-energy. 
Let now na-{k) be discontinuous at k = k^,. With 

^fe.;a:=n.(fc;)-n<,(fe+), (25) 

where k^ is the nearest of the two points k~ and fc^ to 
Sp-cr and (by assumption) Zk^.^r > 0, from Eq. (24) we 
readily obtain 

+ {l--Akt)){el^.,^-el,J. (26) 

For metallic states, where and ef_ are up to in- 

finitesimal corrections equal to Ef for fc^ = kf.„ G S^-a 
[9,10], the LHS of Eq. (26) is vanishing for fc^ e S^-a so 
that, by the arguments presented in [9], Eq. (26) estab- 
lishes the continuity of e^.^ and e^.^ at fc = fc^^ G S^-^a 
(see also § 4.1). 

For fc* at a finite distance from >Sf;<7, the LHS of 
Eq. (26) is strictly positive (see footnote 10), from which 
and from Eq. (26) it follows that at least one of the two 
functions s^.^ and s^.^ must be discontinuous at A; = fc*. 
For fc* sufficiently close to Sp-a, on account of the afore- 
mentioned continuity of e^.^ and e^.^ at all points of 
Sp-cr and on account of the fact that these functions are 
required to achieve the value Sp at the latter points, it 
is most likely that, for a possible discontinuity of e^.^ at 
k = ki,. one has ef^ < ef, and that, for a possible 
discontinuity of at fe = k^,, > e^, . Raising 

the status of these observations to facts, an inspection of 
the signs of the terms on both sides of Eq. (26) reveals 
that the possibility of a continuous e^.^ and discontin- 
uous £^.^ at fe = fe* is ruled out; in contrast, a similar 
inspection reveals that a continuous e^.^ and discontin- 
uous £^.^ at fc = fe* is potentially feasible. On the same 
grounds, it is in principle possible that both e^.^ and e^.^ 
are discontinuous at fc = fe* although, in such an event, 
^k a required to undergo a larger amount of disconti- 
nuity (as implied by Eq. (26)) than is necessary for the 
case in which s^.„ is continuous at fe = fe*. 



It is interesting to point out that the energy disper- 
sion as measured by Zhou et al. [3] corresponding to 
(La2-xSr^)Cu04 in the extreme underdoped limit, where 
X = 0.03, appears to exhibit a finite amoimt of discon- 
tinuity at fe = fe*. On identifying the measured energy 
dispersion with s^.^, it is observed that it indeed satis- 
fies the relationship s"^-,^ < ^^+.^ suggested above on 
general grounds. 

3.5. On the points where e^.^ = e^.^ (and e^.^ = £^.^) 

In our above considerations we have been able to de- 
duce a number of properties associated with e = Sp and 
fe = fepicr G Sp-^a- by specifically relying on the property of 
continuity of e^.^ at fe = kp-^ G <5f;<t (see § 3.4). Here we 
explicitly assume e^.^ to be continuous also at fc ^ Sp-g-, 
specifically at fe = fe* where r\cr{k) is discontinuous; fol- 
lowing the analysis in § 3.4, s^.^ must therefore by neces- 
sity be discontinuous at fe = fe* . Let the latter fe* be fur- 
ther located inside the Fermi sea of the GS of the system 
under consideration; for the ensuing discussions it is not 
necessary that this GS be metallic. The general expres- 
sion for e^.^ can be cast into the form Af/Ag = h{xo) 
(see Eq. (19) above), with xq here identified with fe*, 
/ with the numerator of the expression on the RHS of 
Eq. (1) and g with the denominator (and thus h with 
^fc cr)- From this and Eqs. (8) and (9) it follows that A/ 
is to be identified with Zk^-crS^^ a ^^"^ with Zk^.^a-, so 
that by Eq. (19) we have the following exact result: 

<;. = ^fe*;.' (27) 

where is the exact solution of Eq. (3) and, as men- 
tioned above, fe* a point of discontinuity of no^(fc) at 
which, moreover, e^.^ is explicitly assumed to be con- 
tinuous. Equation (27) is to be contrasted with the more 
specific exact result in Eq. (4) above. In this connec- 
tion we should point out that the validity of Eq. (4) 
does not depend on whether n<^(fe) is discontinuous at 
fe = feF:cr S Sp;a-- This suggcsts that the validity of 
Eq. (27) may similarly not depend on the assumed dis- 
continuity (as opposed to the more general condition of 
singularity) of no-(fe) at the fe = fe* under consideration; 
at the time of writing this paper, wc arc not in a posi- 
tion to make a definite statement on this aspect. In this 
connection note that in general na-(k) is singular (not 
necessarily discontinuous) for all fe G Sp-a [9,10]. 

It is interesting to note that in cases where e^.^ is, 
similar to no^(fe), discontinuous at fe = fe*, application 
of Eq. (22) followed by some trivial algebra yields the 
following exact expression 



With reference to our remark in footnote 9, we point out 
that Eq. (27) is specific to cases in which = 1. 
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= ^(n.(fe;)e<-^^ - n.(fe:)4 J, (28) 

from which one recovers the result in Eq. (27) for e^-.^ 



approaching 



. In cases where na-{k^)/Zk^-a- ~ 1> 
implying thus r\^{k^) « 0, from Eq. (28) it follows that 

^k^:a ~ ^k-■,<T^ 

In the light of the results in Eqs. (4) and (27), it is 
important to point out that, for Fermi- liquid metallic 
states of the single-band Hubbard Hamiltonian (if such 
states at all exist), we have shown [9] that the gradi- 
ent of the quasi-particle energy dispersion S/,.^ (or the 
Fermi velocity times h) for k inside the Fermi sea and in 
the close neighbourhood of iSF:cr differs from '^k£k-a 
a multiplicative constant. A, which we have shown [9] to 
approach unity ^'^ on reducing the on-site interaction en- 
ergy U . It follows that in general the validity of (27) docs 
not extend beyond the points k where n(j(fc) is discon- 
tinuous; in other words, in general the asymptotic series 
expansions of e^.^ and e^.^ centred around a fc, say fe^,, 
at which na{k) is discontinuous (so that e = sat- 
isfies Eq. (3) at fe = fc^; see § 3.1), only agree to the 
leading order and deviate at higher orders in (fc — k-^): 
analogously for £^.^ and £^.^ in the event that e^.^ is 
continuous at fe = fe* (see however the last but one para- 
graph in § 3.4). This aspect becomes clearly evident by 
considering the two-body interaction potential to be the 
long-range Coulomb potential. In this case, Vfce^.^. is 
logarithmically divergent for k approaching the points of 
discontinuity of r\a{k) (see § 4; see also [10]), a property 
that may or may not be shared by V^eJ.^; for instance, 
for Fermi-liquid metallic states, VfeeJ.^ is by definition 
[18] bounded, satisfying 



^k£k„ 



for k^-^fj G Si 



(29) 



The fact that the two important fe points that feature 
in the experimental observations as reported in [3-8], 
namely the nodal Fermi point fcF;CT and fc+, are rela- 
tively very close to one another (typically, Wk^-a — fe*|| 



In [9] for isotropic Fermi-liquid metallic states of fermions 
of bare mass rUe interacting through a short-range poten- 
tial we have obtained A = (m*/mc)nCT(fcj^o.)/[2no-(fe^o-) — 1], 
where m% stands for the ronormalizcd (or effective) mass. 
For the sole purpose of illustration, making the assumption 
that r\a{k^.„) (1 + Zk^.^)/2 and employing the data con- 
cerning ^fep.„ and me/m*^, corresponding to the Coulomh- 
interacting homogeneous electron-gas system, as presented in 
respectively Tables 5.6 and 5.7 of the book by Mahan [23], 
from A « mj(l + Zk^.^)/{2meZk^.^) for (rs;A) we obtain: 
(0;1) [exact], (1;1.12), (2; 1.16), (3; 1.19) and (4; 1.22). Here 
Ts stands for the dimensionless Wigner-Seitz density parame- 
ter, with rs — > corresponding to the uncorrelated limit. 



amounts to approximately 5% of [|feF:crll), implies, follow- 
ing Eqs. (4) and (27), that the fundamental distinction 
between e~^.^ and e^.^ must be of minor observational 
consequence for fe inside the interval between fe^, and 
feF;CT along the diagonal direction of the IBZ compared 
with outside. In other words, in the light of Eqs. (4) and 
(27), it is expected that, for fe inside the latter interval 
(between fe,^ and feF;<7)) ^k-a ^^"^ ^k a should to a good ap- 
proximation coincide. In the specific case at hand, where 
one of our main conclusions drawn from the experimental 
observations in [3-8] , specifically those by Zhou et al. [3] , 
is that the underlying n^{k) is continuous at the nodal 
points of the Fermi surfaces of the compounds studied 
(see § 4.1), the above-mentioned differences between e^.^ 
and £^.^ should be minimal for fe in the close neighbour- 
hoods of the nodal Fermi points. 

Concerning the 'universality' in the nodal Fermi ve- 
locities as observed by Zhou et al. [3], since in all cases 
(whether different cuprate compounds are concerned or 
a specific compound at different levels of hole doping) 
the corresponding Fermi energies are chosen as the ori- 
gin of energy, it follows that the above-mentioned uni- 
versality of the nodal Fermi velocities is indicative of 
the universality of the energy difference Ef — £*, where 
£i,:=e^ = This in turn is suggestive of the pos- 

sibility that the root cause of the discontinuity in r\cr{k) 
at fe = fe* should be external to the electronic degrees of 
freedom, such as longitudinal optical phonons which have 
been suggested [5,8] as inducing the 'kink', or 'break', in 
the observed single-particle energy dispersions in the in- 
vestigated cuprates. For an extensive discussion of this 
aspect see § 5 where we argue that, whereas phonons may 
be vital in bringing about the observed 'kinks' in the en- 
ergy dispersions, they cannot be the immediate cause but 
possibly a secondary determinant through their modifi- 
cation of the electronic states in such a way that n^{k) 
is rendered discontinuous at fe = fe* inside the under- 
lying Fermi seas and that Ef — £* is close to 50 to 70 
meV. It is important to point out that, with the excep- 
tion of rigid parabolic bands, all rigid tight-binding bands 
give rise to relatively strong variation of the bare mass 
of fermions as a function of doping; such variation has 
undoubtedly consequences for the corresponding dressed 
mass and thus the interacting Fermi velocity which, in 
such cases as considered here, can be artificial and un- 
desirable. The experimental observations by Zhou et al. 
[3] are thus suggestive of the possibility that rigid tight- 
binding bands may not be appropriate (specifically) while 
dealing with the cuprate compounds. The unsuitability 
of rigid tight-binding bands (as well as rigid on-site inter- 
action energies) in general is the essential argument un- 
derlying the dynamic Hubbard model proposed by Hirsch 
[24] (see also [25] and the references therein) and their 
shortcomings specifically in applications concerning the 
cuprate compounds has been explicitly emphasized in [9] . 
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For completeness, in cases in which e^.^ is continuous 
at fc = fc* where r\cr{k) is discontinuous (so that, by the 
observations in § 3.4, e^.^ must necessarily be discontin- 
; sec however the last but one paragraph 
in § 3.4), from Eq. (23) and the application of the result 
in Eq. (19) we have (see Eqs. (13) and (27)) 

By relaxing the condition of the continuity of e^.^ at 
fe = fe*, along the same lines as those leading to Eq. (28) 
we obtain 



(31) 



from which Eq. (30) is recovered for approaching 

e^+^. In cases where x\a(k~) / Zk^,;a ~ 1, implying thus 

n<,(fe+) « 0, from Eq. (31) it follows that et « e> . 

* ' k* 

The results in Eqs. (27) and (30), which cannot si- 
multaneously hold (see § 3.4), gain special significance 
by realizing the fact that e^.^ and e^.^ correspond to a 
single-particle spectral function Aa{k;e) introduced in 
[9,10], defined for all k as follows 

A.{k; e):=n{r^^{k) S{e - £< J + [1 - n<,(fe)] S{e - e>„)}. 

(32) 

In [9,10] a number of salient properties of Aa{k;s) have 

been explicitly demonstrated. For instance, on replac- 
ing Acr{k;e) by Aa{k;e) in the expression for the total 
energy of the GS of the interacting Hamiltonian H, one 
obtains exactly the same value, that is the exact GS to- 
tal energy. Although in [9,10] Ac{k;e) was associated 
with some fictitious particles, having some direct signif- 
icance to the properties of the single-particle excitations 
in the neighbourhood of Sf-o-^ Eqs. (27) and (30) unequiv- 
ocally show that Aa{k; e) has special significance also for 
all k outside Sp-^ where r\cr(k) is discontinuous (as we 
have indicated earlier, this significance may also extend 
to all points outside Sp-^ where rt„{k) is merely singular 
and not specifically discontinuous). With this in mind, 
Eq. (32) shows that, at such points as A; = A;,^ where ncr{k) 
undergoes a discontinuous change, this is accompanied 
by an equally discontinuous redistribution of the spec- 
tral weights of the peaks in the single-particle spectral 
function centred at e = e^.^ and e = £^.^ (see Eqs. (27), 
(28), (30) and (31)). Concerning the general subject of 
the spectral weight redistribution in interacting systems 
as described by the Hubbard Hamiltonian, we refer the 
reader to [26], and for related considerations concerning 
specific strongly correlated systems to [27-29]. We em- 
phasize once more that Eqs. (27) and (30) are local (i.e. 
they apply to isolated points, such as k^, is, and not to 

open domains of the k space) so that they do not imply 

< 

equality of 'VkS^.^ with VfcE^.^ at fe = fe*. 



3.6. On the number of solutions of the quasiparticle equation 
at fc = fe^CT corresponding to well-defined quasiparticles 

In § 3.4 we have presented details which rigorously 
demonstrate the continuity of e^.^ for fe in the neigh- 
bourhood of Sp-a (see also § 4) . A corollary to this result 
in conjunction with Eq. (19) is that, in cases where the 



(e) corresponding to fc* 
ishing, one must have 



fcp.^ is not identically van- 



Ef = 



(33) 



This result establishes that, for fe^, = k^-at Eq. (3) can- 
not possess more than one solution corresponding to a 
well-defined quasi-particle; that is to say, in cases where 
ncr{k) is discontinuous at fe = fe* and fe* e >SF;cr, the ex- 
pression for S|]r(e) in Eq. (9) is the most general expres- 
sion. -'^^ The proof of this statement is as follows. Let 
S~(e) = hZi5{e — ei) -f hZ25{e — 82), where £1,62 < M 
and Zi+Z^ < 1 (the following arguments can be trivially 
generalized for any number of similar terms contributing 
to S-(£)). From Eq. (33) it follows that 

_ £1^1 -|- £2-^2 

Z^ + Z2 ' 

which, through the fact that in the case under consider- 
ation £1 = £f (or £2 = £f), implies that Z2 = {Zi = 0). 
This proves our above assertion. In other words, under 
the above-mentioned conditions, Eq. (3) has at most one 
solution at k~ (and similarly for fe^) corresponding to a 
well-defined quasiparticle. 

3.7. On the positivity and the possible discontinuities of 

fe;cr 

The continuity of £^.^ at a point, say fc*. implies that, 
for A£^.^ to be continuous at this point, it is necessary 
and suHicient that the first term in the second expression 
on the RHS of Eq. (7) be continuous at fc = fc*. The 
continuity of the latter term at fe = fe* in cases where 
n„{k) is discontinuous at fe = fe*, following the results in 
Eqs. (19) and (27), imphes that, for Ae^.^ to be contin- 
uous at fe = fe*, one must have 



-(^fc,;.) =0. (34) 



This result is in contradiction with the fact that, for inter- 
acting systems, A£^.^ > (see § 2). It thus follows that 
our assumption with regard to the continuity of both e^.^ 



With reference to our remark in footnote 9, this impUes 
that, for fc* € 5f;<t, < 1. 



10 



and Aff^.^ at fc = fc^,, where no-(fc) is discontinuous, must 
be false. We thus arrive at the general conclusion that, at 
points (which may or may not be points of Sp-a) where 
nr^(k) is discontinuous, e^.^ and Ae^.^ cannot both be where 



continuous; at k 

-< 



where s^.^ is continuous 



F;cr G "^i 

(sec § 3.4), Ae^.^ must be discontinuous in cases where 
n<j(fe) is discontinuous. 



§ 4. Detailed considerations 

4.1. The interacting Hamiltonian and some basic details 

We consider the following interacting Hamiltonian: 

H = Ho+gHi, (35) 

where 

Ho = ^ekcl.^Ck;a, (36) 



k.a 



Hi 



-y 

9.0. ^ 



^k+q;a'-'p-q;(7''-P;° 



o'Ck-a- (37) 



In Eq. (35), g (which has the dimension of energy) is the 

coupling constant of interaction, Sk is a spin-degenerate 
single-particle energy dispersion (which may or may not 
be isotropic), c^.^ and Ck-a are the canonical creation 
and annihilation operators respectively for fermions with 
spin index a, gw(||q||) = t^dlqll) is the Fourier trans- 
form of the two-body interaction potential v{r — r') (as- 
sumed to be isotropic), and Vl = L"^ is the volume of the 
(macroscopic) d-dimensional hypercubic box occupied by 
the system. The wave-vector sums in Eqs. (36) and (37) 
are over a regular lattice (the underlying lattice constant 
being equal to 21: /L) covering in principle the entire re- 
ciprocal space. On effecting the substitutions 



n 



(38) 



and restricting the wave-vector sums to iVt wave vec- 
tors uniformly distributed over the IBZ associated with 
the Bravais lattice spanned by = 1, . . . ,7Vl}, the 

Hamiltonian in Eq. (35) transforms into the conventional 
single-band Hubbard Hamiltonian Ti. corresponding to 
the on-site interaction energy U; in cases where k + q 
and p — q on the RHS of Eq. (37) he outside the IBZ, 
these vectors are to be identified with the vectors inside 
the IBZ obtained from k + q and p — q through Umk- 
lapp processes. On relaxing the replacements in Eq. (38), 
while maintaining the above-mentioned restrictions con- 
cerning the wave vectors, the Hamiltonian in Eq. (35) 
coincides with an 'extended' Hubbard Hamiltonian. 

For the TV-particle uniform GS |4'jv;o) of the Hamilto- 
nian in Eq. (35) we have [9,10] 



ek + g^k;a, Vfe, (39) 



n<7(A;)' 



(40) 



in which 



(^k;a - 77 EE 



a' p',q' 
At At 



X (*JV;o|c;.^Cj^,+q,.^,Cfe+q/;aCp';a' |*JV;0). 



(41) 



We should emphasize that, although e^.^ < fj,, nonethe- 
less is defined over the entire available k space; 
this we have made explicit in Eq. (39) through Vfc. We 
note in passing that, similar to e^.^ and following the 
same techniques as employed in [9,10,19], for systems 
corresponding to bounded £fc, one can readily express 
As^.^{E 00), as defined in Eq. (7) above, in terms of 
GS correlation functions. 

Along the same lines as used in [9,10] it can be shown 
that, for the set <S of fe points for which sf coincides 

fe;o- 

with ffp, the function 



k:a- 



^k:a 



-(Ef 



(42) 



approaches zero for ||fc — fc|[ ^0. This observation is 
significant in that it shows that, for k € S^-^ where 
e^.^ = £f and n„{k) is generically (but not necessarily) 
discontinuous, the energy dispersion e^.^ is continuous in 
the neighbourhood of <Sf;<7 (see also § 3.4). With refer- 
ence to Eq. (19), it follows that, in cases where n^{k) is 
discontinuous at Sp-cr , Pk-a- n^ust undergo a discontinuity 
commensurate with that in n^{k) at S^-cr (and indeed at 
any set S introduced above) for which e^.^ attains the 
value Ep (see § 3.4); in appendix A this aspect is made 
explicit within an approximate framework. In what fol- 
lows, unless we explicitly indicate otherwise, we assume 



'fc;<T 



to be continuous for all k. 



The experimental ARPES data in [3-8] indicate that 
the measured e^.^ (i.e. e"^.^ from our present standpoint) 
behave linearly as function of fe — fep^o- for a finite range 
of fe away from kp-a', focusing on the data in [3], this 
is the case for ||fc — fcF;o-|| varying between zero and ap- 
proximately 5% of ||feF;o-|| for feF:cr a nodal Fermi wave 
vector along the (0, 0) — (tt, tt) direction of the underly- 
ing square-shaped IBZ. With Sk being linear in fe — kp-a 



To tliis end, ouc has to use Eq. (18) in [10], which holds 
for all k, and the fact that and E™(fc), the Hartree-Fock 
self-energy, are continuous for all fe. 
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for such a relatively small interval of k vectors (exclud- 
ing the neighbourhoods of the saddle points of £fe), we 
interpret the experimental observations as implying that 



Cfc;(T ^ b(fcp.^) • (fc - fcp;„) for fc-^fcF;^, 



(43) 



where fepio- is a nodal Fermi wave vector and where we 
assume k to lie inside the Fermi sea. From Eqs. (39), 
(42) and (43) it follows that (see Eq. (64) below) 



With reference to our earlier discussions in § 3.5, we may 
express the Fermi velocity w^,- as encountered within 
the framework of the Landau Fermi-liquid theory as 



(45) 



where A is a constant of the order of unity (see footnote 
13). Experiments by Zhou et al. [3] indicate that vf-_ , 

or v. - very weakly depends on the doping level. 

The observation in Eq. (43) is of far-reaching conse- 
quence. This follows from a combination of two facts. 
Firstly, on k approaching iSpj^ from inside or outside 
the Fermi sea, the following equation has to be satisfied 
[9,10]: 



where 



A<,(fc):= 



V„{k)=K{k), 



n<7(fe) 



■ 1 /I \ ' ro-(K):— — > 
l-n^(fc) efe;a-M 



(46) 



(47) 



For £^.^, defined in Eq. (23) above, we have the following 
equivalent expression [10] 



'fe;cr 



= £fc + 



nSf(fc)-g/3<^ 
1 - n^(fe) ' 



(48) 



where, as mentioned earlier, Y^ili) stands for the 
Hartree-Fock self-energy. It can be shown [10] that, for 
fe kF-a & >5f;o-, one has 



-fe;cr 



/U + a{kF;a) ■ (fc - fepia) + g Ck;a, 



(49) 



£fc;<7 ~ M + a{kF-a) ■ [k - kp-a) - gK{k) {Ck;a- - Vk:a), 

(50) 



where 



a(fep;<,):= VfeEfelfe . 



(51) 



_ nSf(fc)-(£p-£fc,^J 
gnff(fc) 

It can further be shown that [10] 

??fe;(7 ~ for A; ^ fepj^ G iSpj^. 



(52) 



(53) 



Following this and Eq. (50), one arrives at our earlier 
statement concerning e^.^ ^ /i for k — > fepiff G '^f:^; and 
similarly for £^.^, as is apparent from Eqs. (49) and (42) 
(see § 3.4). 

Let T) denote the open domain of wave vectors over 
which nc{k) is continuous. For k & T>, applying Vfc on 
both sides of Eq. (24), we have 

+n„{k)Vke<^ + (1 - n<,(fc)) Vfce>., 

= Vkek + nVkJ:f{k), keV. (54) 

For k in an infinitesimal neighbourhood of S^-a where 
^fc-cr ^'^'i ^fc-<T t^^^ t^*^ value n (up to infinitesimal correc- 
tions), from Eq. (54) we observe that, for a sufficiently 
smooth Efc, a possible singular behaviour in Vfe£^.^ and 
Vfe£^.^ is determined by that in Vk^™{k) on the RHS 
of Eq. (54) . This aspect is explicitly reflected in the 
defining expression for rjk;a in Eq. (52) which through 
Eq. (50) determines the behaviour of £^.^ for A; — > k^-a G 
S-p;ry In other regions of the fc space where e^.^ > s^.^, a 
possible singular behaviour in Vfcl]^''(fe) in some neigh- 
bourhood can in principle entirely or partly be accounted 
for by a similar behaviour in Vfencr(fc). The above 
observations will prove useful in our subsequent consid- 
erations below. 



^® Let n{kp;cr) denote the unit vector normal to <Sf;ct at k^-a £ 
<Sf;ct pointing to the exterior of the Fermi sea FSo-- Since, for fe 
in a close neighbourhood of feF;CT, one has n(feF;CT) '^k^k a < ^ 
for fe G FS<,,FS^ and n(feF;^) • V^e^^^ < for fe G FS<,,FS<,, 
one observes that, in cases where VfcE™(fc) is divergent 
for fe k^-a G <Sf;ct, some useful information can be im- 
mediately deduced from Eq. (54). For instance, for cases 
in which the interaction potential is as long-ranged as the 
Coulomb potential and na{k) is discontinuous at fc = feFio-, 
one can explicitly show that n(feF;o-) • VfcE5^''(fc) — > -|-oo as 
fe — » feF;(T so that, from Eq. (54), one directly deduces that, 
firstly, n(feF;a) • VfcE^.^ ±oo, h{kY:,<j) ■ V^e^.^ =Fcx), 

for fe G FSo-,FSct and, secondly, the balance between these 
diverging contributions must be such that the resultant func- 
tion on the LHS of Eq. (54) approaches -|-oo as fe — > feFicr, 
irrespective of whether fc G FSo- or fc G FSo-. 

In this connection, we point out that according to 
Belyakov [30] (see also [31]) n(fcF;a) • VfenCT(fc) logarithmically 
diverges as fc — » feF;^ for nCT(fe) pertaining to the isotropic GS 
of fermions interacting through a short-range potential. For 
some relevant details see § 6 in [18]. 
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For the interaction potential ?;(r — r') identified with 
the long-range Coulomb potential (or a potential as long 
ranged as this), it can be shown that, so long as Zk^.^ > 0, 
in consequence of E"''(fc) on the RHS of Eq. (52), for 
k — > k~.^ (i.e. for fe approaching kF-a- from inside the 
Fermi sea), one has 

\7Jk,a\ - ||c(fc-J||| In life - fep;,||| life - fep;,||, (55) 

where c(fej;Z^) stands for some constant vector. The com- 
bination of Eqs. (43) and (55) implies that, to leading 
order in (fe — kp-„) for fe — »• k~.^, one has according to 
Eq. (50) 

£fe;a ~ M + g (fe) %;a , (56) 

and consequently 

r^(fe)~0 for k^k-.^ =^ n,(fe-,) = 0. (57) 

Since n„(k^.^) = corresponds to the peculiar condition 
Ac{k^.^;£) = 0, for all e G (— oo,/x] (see Eq. (2) above), 
we discard the condition in Eq. (55) for being unphysical 
in the present context where Eq. (43) is held as valid 
on account of experimental observations. It thus fol- 
lows that, for the case of Coulomb-interacting particles, 
Zfcj,.^ cannot be positive, whereby Eq. (55) is rendered 
obsolete. Assuming Zk^.^ to be positive but relatively 
small, it is not inconceivable that, for fe sufficiently close 
to fepiCT, \Ck;a\ behaves similarly to \r]k-^a\ as presented in 
Eq. (55) which experiments by Zhou et ah [3] may not 
have resolved. Consequently, at this stage the possibil- 
ity of Eq. (55), and thus of Zk^.^ > 0, cannot be un- 
equivocally rejected. We point out that here we are rely- 
ing on our standpoint that experimentally one measures 
£^.^ rather than the 'solution' to Eq. (3). For complete- 
ness, in [10] (§ 6.1.3) the consequences of \Ck;a\ satisfying 
an asymptotic relation similar to that in Eq. (55) have 
been investigated in some detail. In [10] (§ 6.1.3) it has 
been emphasized that, for fermions interacting through 
the long-range Coulomb potential, Fermi-liquid metal- 
lic states belong to the category of states for which, as 
fe — »• fep.^ e Sf-a, rjk-cr satisfies Eq. (55) and (k-a fulfils a 
functionally similar expression. 

We conclude that the assumption of the validity of 
Eq. (43) can be compatible with a picture of fermions in- 
teracting through the long-range Coulomb potential (or 
a potential as long ranged as this) only if Zk^.„. = 0. This 
condition is satisfied (not exclusively) within the frame- 
work where the normal metallic states of the copper-oxide 
based superconductors along the nodal directions of the 
IBZ are of a marginal Fermi-liquid type [32] . This view- 
point enjoys direct experimental support through the 
angle-resolved photoemission experiments by Valla et al. 
[4] along the nodal directions of the IBZ concerning the 
optimally-doped compound Bi2Sr2CaCu208+5. 



The experimental results by Zhou ct al. [3] further 
show that, at fe = fe*, with ||fe* — fepicrll approximately 
5% of ||feF;(7||, fffc-CT is singular, the singularity appearing 
to be a discontinuity in the slope of e^.^ in the overdoped 
regime and some stronger singularity in the underdoped 
regime. The strongest singularity that can be expected 
from e^.^ is a discontinuity (see § 3.4), followed by (in the 
order of significance), for two-body interaction potentials 
as long ranged as the Coulomb potential, a logarithmic 
divergence in Vfec^.^ for fe approaching the points of dis- 
continuity of r\^{k). We believe that the observations by 
Zhou et al. [3] are highly supportive of the viewpoint that, 
at least in the underdoped regime, the long-range of the 
Coulomb potential can be of considerable significance in 
determining the (unconventional) physical properties of 
these compounds. For clarity, for the hole doping frac- 
tion X = 0.03 concerning (La2-a;Sra,)Cu04 (LSCO), one 
observes a behaviour in the measured energy dispersion 
[3] which is reminiscent of a discontinuity superimposed 
by a contribution whose gradient is logarithmically di- 
vergent at fe = fe*. 

4.2. On the relationship between discontinuities in r\a{k) 
and VfeE^.^ 

Here we demonstrate that singularity (not necessarily 
discontinuity) of na-{k) at a point is directly reflected in 
the behaviour of VfeC^.^ at the same point. We shall 
explicitly consider the case where ncr(fe) is discontinuous 
in the interior of FScr\>SF:(T; concerning the behaviour of 
Vfce^.^ for fe fejfjo- e »Sf;(7, this is implicit in our earlier 
work [9,10] which for completeness we recapitulate in the 
following. In this connection, note that r\^{k) is always 
singular (not necessarily discontinuous) on Sp-a- [9jl0]. 

For fe e Si^-a, £k cr and e^.^ attain the value (up to 
infinitesimal corrections) so that the exact inequalities 
^fe-<7 ^ A* and e^.^ > fj, [9,10] imply that, for Sk a smooth 
function of fe in the neighbourhood of S^-a, it is required 
that Vfe^fc:^ undergo at least a discontinuous change 
for fe transposed through iSpjo- from infinitesimally inside 
to infinitesimally outside the Fermi sea. For instance, for 
a specific class of uniform metallic GSs of the single-band 
Hubbard Hamiltonian (to which class the Fermi-liquid 
GSs belong) we have earlier shown [9] (see Eqs. (105) and 
(106) herein) that, with n(k-i,-a-) denoting the outward 
unit vector (pointing from inside to outside the Fermi 
sea) normal to Sp-^ at fe = k^-^ one has 



This aspect is related to the restriction < 7 < 1 for 
the parameter 7 introduced in [9] ; 7 > 1 is excluded for cases 
where > 0. 
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in which 

b- = n{K.,a) ■ Vk^k;a\k=k-„ . (59) 
and (see Eq. (51) above) 

acr:=n(fcp;o-) • a(feF;o-) = n{kF;a) ■ {Jiv'i^l)- (60) 

In [9] it has been shown that the stability of the GS under 
consideration requires that 



b- > 



UK- 



>0, 



where 



1 — n 



> 1. 



(61) 



(62) 



From Eqs. (59) and (58), one observes that on k cross- 
ing 5f;ct (from inside to outside the Fermi sea), indeed 
the derivative of £,k\a with respect to k in the direction 
normal to S^-cr undergoes a discontinuous change (this 
amounts to 2b~ + a^/g > a„{2/A.~ + 1) > 0), from 
the positive value 6~ at fe = k~.^ to the negative value 
— {b~ + Qa/g) at k = fcjjj^. Similar behaviour is shared 
by the (,k;a pertaining to systems of fermions interact- 
ing through two-body potentials of arbitrary range. An 
analogous discontinuity takes place in n{kF-a) ■ ^kS^.^, 
on transposing k through S^-a- For Fermi- liquid metallic 
GSs of systems in which the two-body interaction poten- 
tial is of shorter range than the Coulomb potential, the 
functions e^.^ and £^.^ are such that, on transposing k 

from inside to outside the Fermi sea, the continuous and 

< > 

differentiable extension of £^.^ coincides with e^.^ in the 
close neighbourhood of S-p-^a (sec Fig. 3 in [9]). On the 
basis of this fact and a number of available numerical 
results, in [9] we arrived at the conclusion that, in the 
case of the single-band Hubbard Hamiltonian described 
in terms of the nearest-neighbour hopping integral t and 
the on-site interaction energy U, the corresponding e^.^ 
and e^.^ fail to posses the latter property for U/t^8 (see 
footnote 27 below). 

Having described the singular behaviour of e'^.^ in the 
neighbourhood of S^-a, we now set out to demonstrate 
that a discontinuity in n^ik) at A; = fc* € FSct\>Sf;<7 (i.e. 
fe* is strictly located in the interior of the Fermi sea FS^) 
gives rise to the following result: 



where (see Eq. (13) above) 



''k-y-=l^kek-;a 



fe=fef 



(63) 

(64) 
(65) 



We present the proof of the result in Eq. (63) in § 4.3 
below. Note that < Zk,-cr/f^aiK) < 1- ^® For 
two-body interaction potentials of shorter range than the 
Coulomb potential, Vfe-.^ ~ ^fe+.tj' therefore, on account 
of Eq. (63), one has 



n^(fc-) 



v<+ <2t;f+ . (66) 



'kt-.y 



Introducing, in analogy to the problem of electrons cou- 
pled to phonons [23] (see specifically § 6.4 herein), the 
'mass-enhancement factor' A* as follows 



from Eq. (66) we have 



A, « < 1. 

n<T(fe* ) 



(67) 



(68) 



The result 1 < 1 -I- A* < 2 (specifically 1 + A* w 2) con- 
forms with the experimental observations concerning the 
cuprate compounds [8] (see Fig. 4 herein where A' is what 
here we have denoted by A*) (see also [33]). 



For f fc-.^ we have 



It is important to realize that, in contrast with the case 
where a possible discontinuity in n„{k) at k = fep^o- 6 <Sp;ct 
corresponds to Zk^^ > (see Eq. (16)), there is no a pri- 
ori restriction on the sign of Zk^.^ (viewed as the value of 
na{k^) — na{kt) and not as spectral weight) for k* outside 
<Sf;ct. This is appreciated by realizing the fact that all excita- 
tions corresponding to the wave vector fe* outside Svia can, 
by the very definition of Sp-a, only correspond to many-body 
states whose energies are greater than the energy of the GS 
of the system under consideration, whereby a possible nega- 
tive Zk^-a cannot signify an instability of the latter GS. This 
aspect is already reflected in the fact that, for fc* located at a 
finite distance from 5f:ct, the choice of what we denote by k^ 
and fe^ is in principle arbitrary. That in our present consid- 
erations Zk^-a > Oj is related to the fact that we have chosen 
kt to be the closer of the two vectors fc^ and fcj to the nodal 
Fermi point k^-cr and that, in order for s^,^, or e^.^, to attain 
the required value Sp at fc = feF;^, it must (monotonically) 
increase for k transposed from fe* to k^-a- 
^° We point out that the significance of k^ as the subscript 



of V 



lies in the finite diflierence between e. 



(which lies strictly below /i) and 2^+.^ = (which lies 

strictly above fi) for fc* located at some finite distance from 
iSp;<T (see § 4.3) and not in the possibility that VfcEo-(fe; Sk~ a^ 
is a discontinuous function of fe at fe = fe*. In view of the 
latter, we should emphasize that far from dismissing a dis- 
continuity in VfcSCT(fc; Ej,-.^), at fc = fc*, as being a priori 
infeasible, our statement here only reflects the conflnes of our 
considerations in this paper. We hope to return to this sub- 
ject matter in a future publication. 
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(69) 

whose existence depends on the existence of the second 
term in parentheses on the RHS of Eq. (69). In this 
connection it is important to rcahze that a finite positive 
Zk^-cr is not sufficient for the last term in Eq. (69) to be 
bounded. With reference to the convention adopted in 
Eq. (13), we write 



In a similar fashion, 



-T'^k£k;a\k=kt- 



(70) 



(71) 



For Fermi- liquid metallic states where no-(fc) is discontin- 
uous at fc e «Sp;o-, v^.^ and v^.^ are collinear and point in 
the same direction for all k £ S^-a- For the imiform GSs 
of the single-band Hubbard Hamiltonian, in [9] strong 
evidence is presented indicating that for U/t ^ 8, v^.^ 
and V^.^, k £ S^-.a, arc not collinear so that these states 
cannot be Fermi liquids (see also footnote 27 below). 

Assuming .^|| to be non-zero and bounded, from 
Eq. (63) it follows that i/ ■ differs from u ■ 

for any arbitrary vector f outside the plane normal to 
Vj^ -CT- Let n{ki,) denote the outward unit vector normal 
to iSf;<7 at fepicr S ^f;^, wlth kF-a thc point at which the 
extension of the radius vector fc* meets iSf:ct- For fc* 
inside the Fermi sea and sufficiently close to k^-^a-i the 
inner products of the three vectors in Eq. (63) with h{ki,) 
are all non-negative so that on account of Eq. (63) in 
general we have 



< h{k^) ■ v<+,^ < h{k^) ■ v<_ 



(72) 



The discontinuity in v"^,^ at k = k^,, where by assiimp- 
tion ncr(fc) is discontinuous, is seen to show the same 
trend as observed in the experimental results by Zhou et 
al. [3] (and other workers [4-8]). To clarify this aspect, 
we point out that, by Eq. (4), £^.^ and e^.^ coincide at 
fe = kp-a- Following Eq. (27), whose validity depends 
on the discontinuity of na(k) and continuity of e^.^ at 
fe = fe*, e^.^ and e^.^ also coincide at fe = fe*. This situ- 
ation is schematically illustrated in Fig. 2 where the dis- 
tinctions between the two energy dispersions is depicted 
as being minute over the entire interval between fe* and 
fepjCT along a nodal direction of the IBZ. With this pic- 
ture in mind, the conformity of the result in Eq. (72) (and 
indeed Eq. (63)) with the experimental results by Zhou 
et al. [3] (excluding those corresponding to underdoped 
samples) becomes evident. 

A convenient (though not necessarily an exact) repre- 
sentation of Eq. (63) in terms of two angles and ^ is 



obtained as follows. We take the inner product of both 
sides of Eq. (63) times h with n(fe*) = n(feF;CT) and mul- 
tiply both sides of the resulting expression by 



7:=- 



Ife 



F;a I 



(73) 



so as to obtain a relationship between dimensionless 

quantities. Assuming \\k-p-a — fc*||/||feF;cr|| to be small, for 
similarly small ||fe — fe*||/||feF;CT|| one can consider £^.^ as 
being a linear function of fe, passing, along the nodal di- 
rection of the IBZ, through the point (fe, e) = (feF:cr, Ef) at 
angle 9 with respect to the wave- vector axis (see Fig. 2), 
where 



61 = tan 1 (7n(feF;<,) • Vk£kJk=k 



(74) 



Here < 6 < tt/2, following the convention of count- 
ing angles as positive when they imply counterclockwise 
rotation with respect to point (1,1) in Fig. 2. Follow- 
ing Eqs. (4) and (27), to a good approximation f^.^ 
should be equal to v"^,^ in the interval between fe* and the 
nodal Fermi point feF;^; in the case of particles interact- 
ing through the long-range interaction Coulomb poten- 
tial, the agreement between Vj^,^ and v^.^ in the neigh- 
bourhood of fc = feF;cr is the better the smaller the value 
of Zk^,^, ideally when Z^^^^ — (sec § 4.1). For fc sat- 
isfying ||fe|| < ||fc*||, Vj,.^ and v"^.^ depart according to 
Eq. (63). Since this aspect is a direct consequence of 
^fc*;(T > (see Eq. (63)), it follows that in the neigh- 
bourhood of fe = fe* the behaviour of v"^.^ in relation to 

'"k a ^^^y crucially depends on the range of the two-body 
interaction potential. In view of the experimental results 
obtained by Zhou et al. [3], we believe that the influence 
of the long range of interaction on v^.^ is unequivocally 
present for samples in the underdoped regime; for the 
samples in the optimally-doped and overdoped regimes, 
this aspect is not as unequivocal as in the case of the 
underdoped samples. It is reasonable to believe that this 
feature has its root in the magnitude of Z^^-^cr which for 
an increasing level of (hole) doping should be decreas- 
ing, rendering thus the influence of the long range of the 
two-body Coulomb potential at fe = fe* less effective. In 
order to maintain a non- vanishing amount of discontinu- 
ity in ffc.g. at fe = fe*, however, following Eq. (63), the 
decrease in Zk^-^a- should be accompanied by a concomi- 
tant decrease in no-(fe* ) in such a way that the rate of 
decrease in Zk^;a-/'^a-{k~) (this decrease is evident from 
the the experimental results in [3]) for the increase in the 
level of hole doping is not as strong as would be the case 
for a doping independent ncr{k~). 

For an interaction potential of shorter range than the 
Coulomb potential, it is expected that the behaviour 
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of £^.^ is semilinear 



for sufficiently small values of 
||fe — fc^ll in the region < ||fc*||. This semi-linear line 
stands at angle ^ with respect to the linear line depicting 
e^.^ in Fig. 2, for which we have 



= e- tan 1 (7n(fep;^) • Vkst-Jk^kiJ 



e. 



(75) 



Here — vr < </) < 0, following the convention stated above 
in connection with 9 in Eq. (74). 

In Fig. 3 we present results concerning e^.^ = Sk + 
E,£,k:a within an approximate framework in which ^k-.a 
has the functional form of the Hartree-Fock self-energy 
(see appendix A; Eq. (A21)), evaluated in terms of a 
momentum distribution function n^(k) which is contin- 
uous at the nodal Fermi point fepio- and discontinuous at 
fe = fe*. The results in Fig. 3 concern a planar square- 
lattice model (for the relevant parameters sec the cap- 
tion of the figure) and the two solid curves correspond 
to two densities of particles representing the underdoped 
and overdopcd regions of the hole-doped cupratcs. We 
obtain almost identical results as in Fig. 3 by employ- 
ing an isotropic model defined on the continuum and for 
one, two and three spatial dimensions; our observa- 
tions concerning d = 1 are particularly relevant in light 
of the fact that in the superconducting state, or in the ad- 
vanced stages of the pscudogap phase, the experimental 
observations along the nodal direction of the IBZ should 
be influenced by the reduced wave- vector space available 
to gapless excitations. Not surprisingly, our numerical 
calculations for d = 2 reveal that similarly all qualitative 
aspects of the results in Fig. 3 remain intact by modeling 
r\a-{k) in such a way that the discontinuity in r\a-{k) is 
anisotropic, taking its maximum value at fc = fc* along 
the diagonal directions of the IBZ and diminishing for 
directions away from the diagonal directions. 

Considering the fact that, in real materials, charged 
fermions interact through the long-range Coulomb po- 
tential, the experimental results by Zhou et al. [3], as 



That is, it is linear, but its extension does not pass through 
the 'origin' (1, 1) in Fig. 2. As we shall discuss in some length 
in § 5 below, interestingly the experimental energy dispersions 
in Fig. la of [3], when linearly extrapolated from the region 

ll'^ll < ll'^*!! (i-C- at 'high' binding energies) to fe = k^-a, all 
turn out to meet the energy axis at 75 mcV above the Fermi 
energy Ef- 

22 With t;(||q||) = gw(||q||) and g = eV(47reoerao) (for the 
specification of the quantities encountered here, see the cap- 
tion of Fig. 3), we have employed the following w(||q||): For 
d = 3, M;(||q||) = 47rao/|lqf ; for d = 2, w{\\q\\) = 27rao/||g||; 
for d = 1, w(||q||) = ao exp(||qr||/go) Ei(||q||/go), where Ei{z) 
is the exponential-integral function and go > is a constant 
parameter of dimensions of reciprocal metre. 



well as those by other workers [4-8] , should be viewed as 
consisting of a superposition of the results in Figs. 2 and 
3; here we are leaving aside the possibility of a discon- 
tinuity in e^.^ at fc = fc^, which is theoretically feasible 
(see § 3.4) and appears to be realized in the underdoped 
(La2-xSr^)Cu04 [3], specifically at x = 0.03. 

4.3. Explicit calculation of the amount of discontinuity in 
VfeE^.^ at fe = fc*, with fe* € FSct\<Sf;ct 

Here we present the details underlying the derivation 
of the result in Eq. (63) above. For definitcncss, unless 
we explicitly indicate otherwise, throughout this section 
we assume s^.^ to be continuous at fc = fc* where n„{k) 
is considered discontinuous (sec § 3.4). 

The expression for e^.^ as presented in Eq. (39) is not 
suited for the determination of VfeS^.^, on account of the 
fact that Eq. (39) describes e^.^ in terms of a distribution 
(i.e. Ac{k;£); see text following Eq. (8) above) while the 
available algebraic machinery required in our approach is 
efii'ective for functions. To proceed we therefore define 

ImAk)-=\j^^^z^G,{k;z), m = 0,l, (76) 

in which the contour C has been depicted in Fig. 4. For 
definiteness, we consider C to initiate and terminate at 
the point of infinity in the complex z plane, correspond- 
ing to £^ = cx) in Eqs. (5) and (6). It can be readily 
verified that 



Io;a{k) = n^(fe). 



(77) 
(78) 



Making use of these results, from the defining expression 
in Eq. (39) we have 



Vfe4a = 



VfcJl;a(fc) -£fc;^VfcJ0;a(fc) 

n^(fe) 



(79) 



This expression has the generic form of h{x) defined in 
Eq. (17) so that, in determining the amount of possible 
discontinuity in ^k^k a^ ^® shall rely on the result in 
Eq. (22). 

We assume na{k) to be discontinuous at fe = fe*. In 
what follows we consider the case in which fe* lies strictly 
inside the Fermi sea, that is at a non-vanishing distance 
from Sjr-a- On account of the assumed discontinuity of 
n<^(fe) at fe = fe*, following Eq. (27) we can replace s^.^ 



on the RHS of Eq. (79) by e^-.^ = 



23 



The assump- 



tion with regard to the discontinuity of n^^^k) at fe = fe* 



2^ With reference to our remarks in footnote 9, here we are 
implicitly assuming that = 1. 
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implies that dSg. (k^ ; e) /de is bounded in the neighbour- 
hood of e = £^ (similarly for dSo-(fe* ; e)/de in the 
neighbourhood of e = £^ (see Fig. 1)). Consequently, 
for k in the neighbourhood of fe* we can write 



E<,(fe; z) = S<,(fe; ) + (3k;a {z - e 



for 



fe*:<7/ 



(80) 



1 - ?i/?fe»;a" 

In our subsequent considerations we assume ^k0k;cr to 
be bounded for fc in a neighbourhood of fc* (sec Eq. (94) 
below). For ncr{k) discontinuous at fc = fc* and fe in 
a neighbourhood of fe*, we employ the decomposition 
C = Cs U Cr where Cr = C\Cs', the contour Cg consists 
of the union of two semi-circles of infinitesimal radius 
centred at e^^.^ on the real axis of the complex z plane 
(see Fig. 4). T?hus we write 



W(fe) = /^),(fe)+7W,(fe), 



(82) 



where Im](Tik) and Im-^ik) are the contributions due to 
contours Cs and Cr respectively. The decomposition of 
Im-a{k) according to Eq. (82) can be viewed as arising 
from the decomposition of ^cr(fe; £) in terms of its regular 
and singular contributions, as presented in Eq. (8) above. 
Following our considerations in §§ 3.1 and 3.2, for fe at 
which Eq. (3) has no solution in the interval (— oo,/i), 
Im;a{k) IS identically vanishing. 

From Eq. (79), making use of the general result in 
Eq. (22), we have 

na(fe* ) ^ ^ 

-(Vfe/i;<,(fe) - £fe,;<,Vfe/o;<r(fe)) 

-Zk^;.^ket,X^k-}. (83) 



k=k~ 



In arriving at Eq. (83) we have made use of the assumed 
continuity of £^.^ in a neighbourhood of fc = fc*, imply- 
ing, in conjunction with the assumption of discontinuity 
of no.(fe) at fe = fe*, the equality of with (see 
Eq. (27)). 

The subsequent notation will be considerably simpli- 
fied by introducing the following decomposition 



where (3k;a (not to be confused with (3^.^ introduced in 

Eq. (40) above) is a function which is bounded at fc = fc* 

and o{z — ej^.^^) denotes a function which asymptotically where 

is subdominant with respect to {z — for z — > 

In Eq. (80), T,„{k]z) denotes the analytic continuation 
of Yi„{k\£) into the physical Riemann sheet of the com- 
plex z plane, from which Yj^{k]e) is obtained according 
to I]cr(fe;£) = \iTD.^ioY,„{k;£ ±ivi), £</x. It is trivially 
verified that 

^fe*;a = ^^3— . (81) 



^ fe '^fe;<Tlfe=fe~ 



V*'''£< I + 
^k '^k;a\k=kt' 



o-lfe= 



(84) 



(85) 



(86) 



(87) 



In this way, making use of 

na(fe^) n<,(fe+)' 
Eq. (83) can be written in the following equivalent form 

''fe '^fc;crlfc=fe7 

rio- fe^ M- ^ ^ ife=fe» 



Vfc/l;a(fc) -£fc^.„VfeV(fc) 

Zk^,;a^k^ Sk;a\k=k~\ 



\k=k+ 

(88) 



We now proceed with the determination of the terms 
on the RHS of Eq. (88). It is readily verified that (see 
footnote 24) 



Vfe4;^(fe)|fe=feT = Vfen,(fe)|^^^j , 



(89) 



^'^^^('^) u^uT = (^'^"-('^)lfe=fef 



+ (vr4Jfe=fef )n.(fe?). (90) 

Following the Dyson equation, making use of Eqs. (80) 
and (81), we have 



VkG„{k;z) 



^ Ji -^fe*;o- 
fe=fe- {z-£k-;ay 



■^fe-;a(^). 



for z^e^-.^, (91) 



Below (as well as elsewhere in this paper) =F in 
^fc/(fc)lfe=feT denotes left/right gradients of /(fc), defined as 
the limit of Vfe/(fe) for fe approaching fc* from the left/right; 
here left /right is defined by n{ki,) ■ (fc — fc*) > 0. For the def- 
inition of n(fc*) see text preceding Eq. (72) above. 



Note that here Vj^^'-'e^.^, x = r, s, stands for an entire sym- 
bol, that is V^''^ does not on its own denote an independent 
operation. 
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where (see Eq. (69)) 



1 



k=kZ 



(92) 



Thus, on the basis of Eqs. (76) and (91) and the Cauchy 
theorem, we have 



Vfe/^Ufc) 



k=k- 



C, 2TTi (z-efe;-;^)2 



From Eqs. (92) and (80) we obtain 



d 



Zkticr '^kPk;a\k^ 



SO that from Eq. (93) it follows that 

^Zk,-a ^k(3k-a\k=k- 



4:^(fc) 



k=k~ 



(93) 



(94) 



(95) 



Vfc /g(fe) 



T^Zk^-0 



-x{^k--A^k-;a) + Zk.W^k-;. kfikAk^k') ' (9^) 

Combining the results in Eqs. (89) and (90) we have 

(Vfe/£(fc)-£^,^,V./«(fe))|^^^^^ 

= (Vi-^^^jfe^fej) n.(fej), (97) 
while combining Eqs. (95) and (96) we obtain 



(98) 



where we have made use of (see Eqs. (69) and (92)) 

^fe,:<T =^fe-;a(^ = £fc;;a)• 
Evidently (see Eq. (86)), 



(Vfe/g(fe)-e^^^,Vfe4J(fe)) 



fc=fe+ 



= 0. 



(99) 



From Eqs. (88), (97), (98), (99), making use of Eqs. (82) 
and (85) we arrive at 



< I _ nZk,;a 

^ k ^fc;crlfe=fe- — „ /.-X "k^,■,a^ 



(100) 



which is simply the expression in Eq. (63) in disguise; this 
is readily verified by employing the results in Eqs. (85) 
and (86). We observe that a finite Zk^-a gives rise to a 
finite discontinuity in 'Vk£k-„ at fc = fe* (sec Fig. 2). 

An important aspect of the result in Eq. (100) is the 
following. In § 4.1 we have emphasized that, for metallic 
GSs of fcrmions interacting through potentials which are 
as long ranged as the Coulomb potential, Vfc£^.^ loga- 
rithmically diverges for k approaching any region where 
na{k) undergoes a discontinuous change [10]. Since 



v^^V< I 

''fe '^k;(7\k=k~ 



Vfe4Jfe=fe- -Vfee<.Jfe^fe+, (101) 



in follows that any logarithmically divergent contribution 
that may be contributing to Vfe£^.^ for fe — > fej, which 

by necessity has the same sign for k ^ k~ and k ^ k'^ , 
^® cancel. Thus, even though for the above-mentioned 
systems VfeE^.^ is divergent for k ^ kf, for these sys- 
tems, following Eqs. (100) and (101), is bounded. 

Finally, we point out that Eq. (100) is not applicable 
to cases in which fc^^ € <SF;cr. This aspect is rooted in the 
failure of Eq. (85) for k^, G Sp-^ which, in turn originates 
from the fact that, on all points of Sp-cr, £k ,7 acquires its 
maximum value (namely £p) [9,10]. As can be explicitly 
shown, for fc* G S^-^, the LHS of Eq. (85) is equal to 

^k ^k;a\k=kt- 



§ 5. Summary and discussion 

In this paper we have presented a set of complemen- 
tary pieces of evidence indicating that the single-particle 
energy dispersions as deduced through the ARPES are in 
fact e^.^ (Eq. (39)) rather than ek;a (more precisely, s'^.^) 
that formally is a solution of the quasi-particle equation 
(Eq. (3)). We have shown that, in addition to e^.^ co- 
inciding with £^.^ on the Fermi surface S^-a- of metallic 
GSs, the two energy dispersions are equal at any point 
inside the underlying Fermi sea FSo- where the GS mo- 
mentum distribution function r\a{k) undergoes a discon- 
tinuous change and e^.^ is continuous; at the time of 
writing this paper the question remains open whether 
this property is retained when the last-mentioned dis- 
continuity in ncr(fc) is replaced by a more general kind 
of singularity. One of the conclusions that we draw from 
the analysis of the experimental data obtained by Zhou 
et al. [3] as well as by other workers [4 -8] is that (in par- 
ticular), along the nodal directions of the IBZs of the 
cuprate compounds investigated by the latter workers, 



^® See footnote 16. Further, for a simple model, such as 
the electron-gas model, this aspect can be demonstrated by 
explicit calculation. 
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na{k) is discontinuous at a point k^, located in the close 
vicinity of the nodal Fermi point fepio- (||^f;<t — ki,\\ typ- 
ically amounts to some 5% of ||feF;o-||)- Prom this and 
our aforementioned finding with regard to the equality 
of £^.^ and £^.^ at fc* and kF-a we have concluded that 
in the interval along the diagonal direction of the IBZ 
between k^, and k^-a the two energy dispersions e^.^ and 
^k <7 should be in general experimentally indistinguish- 
able (see Fig. 2). For fermions interacting through the 
long-range Coulomb potential, e^.^ and e^.^ have dis- 
tinct behaviours in the neighbourhoods of points where 
n^(fc) is discontinuous (§ 4.1). Since, according to our 
analyses, the n^{k) corresponding to systems studied by 
Zhou et al. [3] must be either continuous or only slightly 
discontinuous at fe = k^-a, we have arrived at the conclu- 
sion that the influence of the long-range of the Coulomb 
interaction potential should be strongest at fe = fe*. The 
experimental data obtained by Zhou et al. [3] indeed 
show that for the underdoped sample (La2-2;Sr,j;)Cu04, 
corresponding to for instance x = 0.03, the measured 
energy dispersion shows behaviour in the vicinity of fe* 
that markedly differs from similar energy dispersions cor- 
responding to optimally-doped and over-doped samples 
of (La2-a;Sra;)Cu04. This difference can be viewed as 
consisting of two separate contributions, namely a flnitc 
discontinuity and a logarithmic divergence in the slope 
of the energy dispersion at fe = fe*. The former is consis- 
tent with the exact result obtained in § 3.4 according to 
which £^.^ can undergo a discontinuous change at points 
inside FSo- where no-(fe) is discontinuous, and the latter 
is feasible only if the two-body interaction potential is 
as long ranged as the Coulomb potential. We point out 
that, in cases where e^.^ is discontinuous at the point of 
discontinuity of no-(fe), that is fe = fe*. neither t^_^ nor 

e^+^ coincides with e^^.^; here ||fe* 1| < l|fe* 1|- However, 
in cases where r\^{k^)/Zk^-^ » 1, and thus ncr{k^) « 0, 
it can be rigorously shown that e^-.^ ~ ^fe*;cr- 

Our point of departure in this paper has been based on 
the consideration that single-particle energy dispersions 
as measured by the ARPES coincide with the average 
value of e with respect to the normalized energy distri- 
bution function A^{k; e)/r\^{k); this average value 
we have denoted by s'^.^. Throughout this paper we 
have presented a host of evidence to support this point 
of view. On the basis of this viewpoint we have pro- 
posed the standard variance Ae^.^ of e with respect 
to the above-mentioned normalized energy distribution 
function as corresponding to the measured width of the 
peak in Aa-(k;s) centred at s'^.^. From the expression 
for Ae^.^ we have inferred that sharp variations in no-(fe) 
arc directly reflected in the behaviour of Ae^.^. We have 
shown for instance that the combination of a discontin- 
uous no-(fe) and continuous e^.^ at fe = fe* (whereby 
^fe* <T ~ ^kt-a) iniplies a discontinuity in Ae^.^ at fe = fe* 



with the change in Ae^.^ being of the same sign as that in 
n,(fe),thatisn,(fe;)<n,(fe+) ^ Ae<_^^ < A£<^.^. 

For (La2-xSra,)Cu04 at a; = 0.063, Zhou et al. [3] ob- 
served a decrease in the electron scattering rate (i.e. the 
width of the MDC) at the same energy as a discontinu- 
ity in the gradient of the energy dispersion is observed 
on transposing fe from 'below' to 'above' fe*, correspond- 
ing to ||fe|| < ||fe*|| and ||fe|| > ||fe*|| respectively. This 
decrease in the scattering rate is in full conformity with 
the sign associated with Zk^;a-—ncrik^) — r\„{k'^) which 
has to be positive in order for the change in the energy 
dispersion as described by the expression deduced in this 
paper (Eqs. (63) and (100)) to conform with the exper- 
imentally observed change (sec footnote 19). A sudden 
change in the spectral width at the energy corresponding 
to the 'break' in the single-particle energy dispersion has 
also been reported in [5,7,8]. 

Since J2k a ^kr^cr{k) amounts to the kinetic energy cor- 
responding to the GS of the Hamiltonian H in Eq. (35), 
one observes that a finite discontinuity in r\a-{k) in the 
interior of the underlying Fermi sea, of the type indi- 
cated above, necessarily results in a higher kinetic energy 
than the kinetic energy corresponding to a comparable 
na-{k) which, however, is free from the latter discontinu- 
ity; in this connection, recall that no-(fe) > and that 
J2k cr '^crik) = N, independent of the the strength of the 
interparticle interaction, so that a discontinuity in na{k) 
over a set of points {fe*} in the interior of the Fermi sea, 
with nCT(fe* ) > n(j(fe^), is necessarily accompanied by a 
'spillage' of (fe) into the region outside the Fermi sea. 
We thus observe that the interpretation as put forward 
in the present paper of the experimental observations in 
[3-8] also provides a natural explanation concerning the 
experimental observations of the excess 'kinetic' energy 
in the normal states of the doped cupratc compounds 
[11,12] in comparison with the 'kinetic' energies of the 
GSs of conventional metals; in arriving at the latter state- 
ment wc have relied on the supposition that the gaplcss 
excitations along the nodal directions of the IBZ (below 
the pertinent superconducting transition temperatures) 
are akin to the gapless excitations of the normal states of 
the investigated cuprate compounds. Experimental ob- 
servations, at temperatures elevated above the pertinent 
superconducting transition temperatures, along various 
directions of the IBZ [5] (concerning change in the spec- 
tral widths see [8]), are in conformity with the above- 
mentioned supposition. 

Although the considerations in this paper provide a 
consistent scenario concerning the main aspects of the 
experimental data in [3-8], a fundamental question re- 
mains to be answered, namely whether a single-band 
model which, moreover, exclusively takes account of the 
electronic degrees of freedom (Eq. (35)), allows for a 
uniform metallic GS whose corresponding single-particle 
spectral function involves a coherent contribution below 
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Sp characteristic of a well-defined quasi-particle excita- 
tion at fc = fc^ located inside the underlying Fermi sea 
(i.e. at some finite distance from Sp-a)- Note that, al- 
though the existence of a well-defined quasiparticle exci- 
tation at fc = fc* is sufficient for rendering na{f^) discon- 
tinuous at fe = fc*, this discontinuity does not necessarily 
imply that the pertinent excitation energy should lie by 
a finite amount below the Fermi energy; conversely, in 
cases where the excitation energy at issue is less than the 
Fermi energy by a non- vanishing amount, it is necessary 
(in fact by definition) that ki, be located away from <Sf;<7 
(see paragraph following Eq. (12) above). Recent Monte 
Carlo results corresponding to projected variational wave 
functions [34,35] concerning the conventional single-band 
Hubbard Hamiltonian (involving only an on-site interac- 
tion between fermions of opposite spin) on a square lat- 
tice, are not supportive of the possibility of a noticeable 
discontinuity in no-(fe) in the interior of the underlying 
Fermi sea. Furthermore, although for the hole-doping 
fraction x approaching zero, the discontinuity of no-(fc) 
at the nodal Fermi point fep-.tr diminishes (according to 
the authors of [34,35], Z^^,^ scales like x for sufficiently 
small x), for x corresponding to optimal doping (here 
a^opt = 0.18) it is relatively too large for the following 
two reasons. Firstly, since < ncr(fc) < 1, the larger 
the value of Zk^,^ , the smaller is the maximal value that 
■^fe,;<7/no-(fe^), or A* (see Eq. (68)), can possibly attain; 
following Eq. (63), a small value for the latter quantity 
corresponds to a small 'kink' in the single-particle en- 
ergy dispersion e^.^ at fe = fe*. Secondly, in principle 
a non- vanishing (and in practice a relatively large) Zk^,.^ 
corresponding to a nodal Fermi point at optimal dop- 
ing does not conform with the experimental observations 
by Valla ct al. [4] according to which the self-energy in 
the neighbourhood of a nodal Fermi wave vector along a 
nodal direction of the IBZ satisfies the scaling behaviour 
characteristic of the sclf-cncrgics of marginal Fermi liq- 
uids. Assuming that, insofar as our present considera- 
tions are concerned, the variational results in [34,35] are 
sufficiently accurate, we have no alternative but to pro- 
nounce the conventional Hubbard Hamiltonian (as op- 
posed to 'extended' Hubbard Hamiltonians and as op- 
posed to Hamiltonians, such as the Hubbard-Holstein 
Hamiltonian (see, e.g., [36]), in which the fermionic field 
is coupled to an external bosonic field) as being inade- 



■^'^ It is interesting to note that the na(k) as reiiorted iu 
[34,35] (along the nodal direction of the IBZ of a square lattice 
and corresponding to a projected variational wave function 
for U/t = 12 and t' = t/4) unequivocally violates the result 
[ricrikr-a) + "a{kt;a)]/'2 = 1/2- FoUowiug the analysis in [9] 
and in full conformity with the conclusion arrived at herein, 
we observe that the underlying metallic GS is not a Fermi 
liquid. 



quate for describing the phenomena observed in [3] and 
[4-8]. This conclusion contradicts the finding by Rande- 
ria et al. [37] whose analysis is based on the very same 
computational results (i.e. those reported in [34,35]) as 
relied upon by us above. 

Concerning the interpretation of the experimental ob- 
servation in [3-8] as presented in [37] , denoting the Fermi 
velocity at high binding energies by Vp^^^ (to be compared 
with what we have in this paper denoted by v^-.^) and 

that at low binding energies by (to be compared with 

t!<+^), Randeria ct al. [37] deduce that v^'^^ = vl°^/Z 

(compare with Eq. (63) above) where Z stands for what 
we have in this paper denoted by Z^^,^, to be strictly 
distinguished from Zk^.a- The expression as deduced in 
[37] for Vp°™ is the standard expression for the Fermi ve- 
locity in (isotropic) Fermi-liquids (sec Eq. (7) in [37]). 
Through stating that "at intermediate to high energies, 
\dT,'/du)\ <C 1", where S' is the real part of self-energy 
as evaluated at "fc = fep", Randeria et al. [37] deduced 
an expression for v^^^^ which up to the numerical factor 
Z (by the latter statement this is identified with unity) 
coincides with that for i>p°^ and hence the above-quoted 
expression relating the two velocities. It is evident that 
the crucial, and in our opinion manifestly unjustified, 
step taken by Randeria et al. [37] in deducing the ex- 
pression Vp^^^ = uJP^/Z is that of effecting the condition 
"\d'S,'/duj\ ^ 1" in the expression concerning fp'^*^, for 
the essential aspect of the experimental observations be- 
ing interpreted concerns a change in energy dispersions 
(or discontinuity in Fermi velocities) over a relatively very 
narrow range of energies (centred around, for instance, 
approximately = £f — 70 meV insofar as the experi- 
ments reported in [3] are concerned) whereby an appeal 
to the relationship "|9S'/(9w| ^ 1" is seen to be entirely 
unwarranted. The only possibility that would justify an 
appeal to the latter relationship would be the disconti- 
nuity of "dYi'/doj" at — with S' reducing to 
a nearly constant function of oj for all uj < io^. By 
the Kramers-Kronig relation the imaginary part E" of 
the self-energy would be almost vanishing for all < w*, 
which is neither realistic nor enjoined by the experimen- 
tal observations under consideration. We therefore con- 
clude that the considerations by Randeria et al. [37] do 
not shed light on the observations as reported in [3-8] . 

It is important to point out that according to a widely, 
but not universally, held view (for a review see [33]; 



For the uniform GSs of the single-band Hubbard Hamil- 
tonian, which are dealt with in [37], lim|e|^oo Sa-(fc; e) = 
fi~^Un^, where ria is the number of fermions with spin index 
(T (the index complementary to a) per site and U is the on- 
site interaction energy. The above-mentioned nearly-constant 
function must therefore be close to the latter limiting value. 
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see also [38]), both the 'kink' in the single-particle en- 
ergy dispersion and the decrease in the spectral width at 
fe = fe* are to be attributed to electron-phonon interac- 
tion [8,39]. Counterarguments, disfavouring phonons 
and attributing the latter phenomena to antiferromag- 
netic spin fluctuations have been put forward [43,44,42] 
(see also the discussions in [6,7] which concern exper- 
iments as well as those in [5], examining experimental 
observations in the light of a number of different theoret- 
ical possibilities). We do not consider electron-phonon 
interaction as being the primary cause for the 'kink' in 
the experimentally measured energy dispersions. This 
becomes evident by linearly extrapolating the data in 
Fig. la of [3] at high binding energies (corresponding 
to energies between approximately 0.2 and 0.1 eV be- 
low the Fermi energy Sp) towards the nodal Fermi point 
fe = kp-a- One observes that at fc = k^-cr the extrapo- 
lated energy dispersions corresponding to all doping con- 
centrations (in total 10 energy dispersions correspond- 
ing to a; G [0.03,0.3]) intersect the energy axis at the 
same point (this with remarkable accuracy) above the 
Fermi energy, namely at approximately 0.075 eV above 
Ef (instead of at, or nearly at, Sp if electron-phonon 
interaction were directly at work; see later); this ex- 
cess energy is almost equal to the binding energy (i.e. 
Ef — £*) at which 'kink' in the energy dispersions has 
been observed. The conventional electron-phonon the- 
ory [23] (sec in particular § 6.4 herein), according to 
which ||ffc+.^||/||t^fc-.<,|| = 1/(1 + A), where A is the cou- 
pling constant of the electron-phonon interaction, or the 
electron-phonon contribution to the mass enhancement 
(for a comprehensive overview see [45]), requires that the 
change in the energy dispersion correspond to a clockwise 
'winding' around Sp (inside an energy window of width 
hojph centred at e^, where Uph the the characteristic 



Concerning the role envisioned for the electron-phonon in- 
teraction in the context of high-temperature superconductiv- 
ity in the cuprate compounds see [40] and for the interplay be- 
tween electron-electron and electron-phonon interactions see 
[41]. 

^° For completeness we mention that we fail to comprehend 
the methodology with the aid of which the authors of [42] de- 
duced from their computational results concerning the singlc- 
particlc spectral function (denoted in [42] by N{k,ij)) a kiuk 
in the single-particle energy dispersion. In our judgement, 
connecting the maxima of a series of N{k,u)) corresponding 
to various k (considered as parameter), as the authors of [42] 
did (see Figs. 4 - 6 in [42], noting that in these ui is measured 
in milielectron volts and k in units of the inverse lattice con- 
stant), cannot provide unequivocal information with regard 
to the dispersion of the underlying single-particle excitation 
energies. 

See, for example. Fig. 26.1 in [46]; see also Fig. 6.16 in [23] 
and note that the electron-phonon self-energy is vanishing at 



phonon frequency) of the electronic dispersion with re- 
spect to that in the absence of phonons; the 'kink' in 
the single-particle energy dispersions in systems of elec- 
trons coupled to phonons is thus seen to correspond to 
a sharp crossover of the phonon-dressed dispersion into 
that of electrons in the absence of electron-phonon in- 
teraction (from which the phonon-dressed energy disper- 
sion becomes rapidly indistinguishable) at some short 
distance outside the above-mentioned energy window; 
linear extrapolation of electron-phonon dressed energy 
dispersions from outside the above-mentioned crossover 
region (and corresponding to \\k\\ < ||fe*ll) must there- 
fore at fe = fcF:cr very nearly intersect the energy axis at 
Ep. Experimentally, this aspect is evident from the high- 
resolution photoemission data concerning the (0001) sur- 
face state of beryllium (also referred to as the F-surface 
state of Be(OOOl)) as presented in Fig. 3 of the work 
by Hengsberger et al. [47] (see also [48]). A comparison 
of the data in Fig. 3 of [49] with those in Fig. 4 of [48] 
leads one to conclude that the above-mentioned clockwise 
'winding' around Sf of the phonon-dressed single-particle 
energy dispersions with respect to their undressed coim- 
terparts is also present in the energy dispersion of the 
(110) surface-state of molybdenum [49]. With reference 
to our above statement with regard to the extrapolation 
of the energy dispersions in Fig. la of [3] corresponding 
to life]] < IJfc*!] towards fe = fcF;cr, it is evident that the 
'kinks' observed in the single-particle energy dispersions 
in the hole-doped cuprate compounds [3-8] cannot be a 
direct consequence of the electron-phonon interaction. 

In order to examine the possibility of an indirect role 
played by phonons in bringing about 'kink' in the energy 
dispersions as observed in [3-8] , it is instructive to recall 
a model proposed by Frohlich [50,51] (see also §§ 1.8 and 
3.16 in [52]). Frohlich [50,51] showed that, for sufficiently 
large value of coupling constant, the interaction of inde- 
pendent electrons with phonons gives rise to an exotic 
electronic GS, which Frohlich showed to be a supercon- 
ducting state. In this state, the Fermi sea of electrons 
consists of two concentric regions, separated by a narrow 
region of unoccupied single-particle states (see Figs. 1 
and 2 in [50]); Frohlich explicitly showed that this GS, 
corresponding to a 'shell' distribution of electrons in the 
fc space, is characterized by the equality of the widths of 
the intervening unoccupied region and the outer occupied 
region of the fe space. For such a state, on transposing 



e = Ef, or using the notation adopted in [23], at u = 0. 

Schafroth [53] has, however, shown that, as the theory 
of Frohlich is based on a finite-order perturbation theory 
concerning electron-phonon interaction (explicitly, a second- 
order theory), it fails to describe the Mcissner effect. For an 
exposition of the work by Schafroth see [52] (appendix A1.8 
herein). 
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k from the interior of the inner Fermi sea to the exterior 
of the outer Fermi sea, ncr{k) undergoes three disconti- 
nuities, from unity to zero, from zero to unity and again 
from unity to zero respectively. Each of these discontinu- 
ities would in turn give rise to discontinuities in e"^.^ (in 
principle; § 3.4) and ^k£k,r, in the manner described in 
this paper (§§ 3.4, 4.2 and 4.3), were it not for the fact 
that these discontinuities all correspond to single-particle 
excitations whose energies coincide with Sp. It follows 
that the Frohlich state is not the appropriate state inso- 
far as the observations in [3-8] are concerned. In spite of 
this, the crucial role that phonons play in bringing about 
the Frohlich state (at least in principle), with the above- 
indicated discontinuities in the associated na{f^), makes 
evident that phonons can in principle be vital in regard 
to the observations reported in [3-8]; in this connection 
it should be noted that the binding energies £f — £* corre- 
sponding to 'kinks' in the measured single-particle energy 
dispersions of the investigated cupratc compounds very 
nearly coincide with the energies of the longitudinal op- 
tic phonons in these [8] (for a detailed discussion of this 
subject sec [33]). 

We point out that the prediction by Frohlich [50] of 
the above-mentioned 'shell' state relies on a second-order 
perturbation expansion of the GS energy of the coupled 
electron-phonon system. Soon after its predication, this 
state was asserted by Wentzel [54] and Kohn and Vachas- 
pati [55] as not feasible through its interception by a 
structural lattice instability at an electron-phonon cou- 
pling constant smaller than that required for the forma- 
tion of the Frohlich state; Wentzel [54] conceded, how- 
ever, that interaction of electrons, not taken account of 
in Frohlich's considerations, would render the Frohlich 
state feasible. However, according to Bardecn [56] (see 
the closing paragraph of section III herein), who has es- 
tablished a modified stability condition for the lattice, 
the Frohlich state is viable. 



Anderson [57] presented the eontents of the Frohlich 
model in a modern setting, indicating the diagrammatic rep- 
resentations of the underlying exchange processes. Subse- 
quently, with reference to a Migdal theorem [58], Anderson 
[57] stated (p. 1346): "But in fact Migdal resolved the prob- 
lem even more completely by showing that E is a much more 
sensitive function of oj than of k — i.e., it is local in space, 
retarded in time, thus depends on co much more sharply — 
so that the large correction to (-Efc)off comes from dTi/du>: ... 
and Z — — (dTj/dij)] is always positive: the Frdhlicli- 
Bajrdeen instability simply does not occur. The complete 
treatment shows that phonon instability and the first sin- 
gularity of Z actually occur precisely at the same coupling 
strength." Two comments arc in order. Firstly, the electronic 
state on which the Migdal theorem at issue is based is that 
of the non-interacting free-fermion model (this is readily veri- 
fied through inspecting [58] in conjunction with [50]); in other 



Concerning the 'universality' (within an experimental 
uncertainty of approximately 20%) of the Fermi veloci- 
ties at low binding energies as reported by Zhou et al. [3] , 
following the analysis in § 3.5, this 'universality' reveals 
two aspects; firstly, the 'universal' energy — where 
e* = = is associated with some bosonic ex- 

citations that are external to the electronic degrees of 
freedom, supporting the possibility that longitudinal op- 
tical phonons, though not directly (see above), can indi- 
rectly be the sought-after bosonic excitations; secondly, 
a rigid, that is doping-independent, tight-binding energy 
dispersion Sk underlying the many-body Hamiltonian in 
Eq. (35), is not capable of reliably describing the cuprate 



words, it is assumed that the scale of the electronic excitation 
energies is determined by a single electronic mass parameter. 
For specifically strongly-correlated electron systems there is 
no a priori reason to believe that the GS of this model would 
at all bo an appropriate starting point for dealing with the 
problem of electron-phonon interaction (this is exactly the 
same problem which in fact hinders a rigorous formulation of 
a theory concerning the superconducting states of the cuprate 
compounds whose low-lying single-particle excitations in the 
normal state are not quasi-particle-likc; see however [59]), so 
that the above-mentioned Frohlich-Bardeen instability can- 
not on its own rule out the existence of the Frohlich state. 
Secondly, it is generally, even though perhaps not universally, 
accepted that non-adiabatic processes are significant in the 
copper-oxide-based high-temperature superconductors, lead- 
ing to the 'failure' of the Migdal theorem referred to above. 
See [60]. For completeness, we have carried out calculations 
[61] on a uniform two-dimensional model in which the GS is 
the Frohlich state. From the first-order self-energy in terms 
of the dynamically screened interaction function we have de- 
duced quasi-particle lifetimes which are generically by one or- 
der of magnitude shorter than the lifetimes of quasi-particles 
in the conventional two-dimensional model at similar densi- 
ties. We have further studied [61] a variety of the properties of 
this model in the superconducting state within the Bardeen- 
Cooper-Schrieffer (BCS) framework, despite the fact that the 
aforementioned short quasi-particle lifetimes undermines the 
use of the BCS formalism. For instance, in the weak-coupling 
limit we obtain that 2A(T — 0)/[A;bTc] is relatively close to 
4 (to be compared with the conventional BCS value of 3.52). 
We point out that a recent [62] non-selfconsistent first-order 
calculation (in which the effective interaction has been de- 
duced from diffusion Monte Carlo results concerning the static 
charge and spin structure factors and random-phase approxi- 
mation Ansatze concerning the dynamical density-density and 
spin-spin susceptibilities) concerning two-dimensional liquid 
''He indicates that on increasing the aerial density of ''He 
atoms the effective mass diverges; a subsequent increase of 
the density suggests what Boronat et al. [62] refer to as being 
indicative of a transition to "anomalous occupation number^' 
and to which we would refer as the Frohlich 'shell' state. In- 
terestingly, according to Boronat et al. [62] "Both spin and 
density fluctuations have profound effects." 
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compounds, specifically those investigated by Zhou ct al. 
[3] ; the bare mass corresponding to a rigid tight-binding 
energy dispersion is relatively strongly dependent on the 
level of doping (sec caption of Fig. 3), and this is most 
likely inherited by the corresponding dressed mass and 
consequently the interacting Fermi velocity, thus contra- 
dicting the above-mentioned experimental observations. 

Although in this paper we have mainly focused on the 
ARPES data along the nodal directions of the IBZs of the 
cuprate compounds as reported in [3- 8] , our approach is 
applicable along any direction of the IBZ. In the light of 
the observations as reported in [5] (see specifically Fig. 3 
in [5]) we conclude that a discontinuity in n^{k) in the 
interior of the Fermi sea may be a generic property of the 
normal states of all cuprate superconductors. With refer- 
ence to our observation in this paper that along the nodal 
directions of the IBZs of the cuprate compounds inves- 
tigated in [3-8], the quasi-particle equation must have 
at least (see footnote 9) two solutions < /i and 

^fe*-cr ^ A* corresponding to true Landau quasi-particles 
(in the sense of Zk^-a- > 0), it is interesting to recall an 
earlier observation presented in [9]; here it has been rig- 
orously shown that, within the framework of the conven- 
tional single-band Hubbard Hamiltonian, for k located in 
the pscudogap regions of the reciprocal space, the quasi- 
particle equation must have (at least) two solutions, one 
strictly below n and one strictly above n, and that these 
solutions do not correspond to quasi-particle excitations 
of the Landau type, but to resonances. It is of interest 
to investigate whether, and if so, how these two distinct 
energies are related to and which correspond 

to Landau quasiparticles. We should emphasize that, 
whereas the above-mentioned resonances correspond to 
wave vectors in the pscudogap regions of the underlying 
putative Fermi surfaces (in [9] it has been shown that 
S^^^XSp^a, which may or may not be empty, constitutes 
the pscudogap region of the wave- vector space), the wave 
vector fe* is located in the interior of Fermi seas, i.e. at 
some finite distance from <Sf;o-- 

We close this section by reiterating that, although we 
have been able to provide a consistent scenario for the ex- 
perimental observations in [3-8] , a number of outstanding 
questions remain. These are 

(i) whether a single-band model is capable of possess- 
ing a uniform metallic GS whose momentum distribution 
function is discontinuous in the interior of the underlying 
Fermi sea (specifically along the nodal directions of the 



Here Sp-l stands for the Fermi surface corresponding to 
+ ia, in which ta has been chosen such that the interior 
of Sp°^ contains the same number of k points as contained in 
the interior of S^^a, namely N^, the total number of fermions 
in the GS with spin index a. 



IBZ) and at the pertinent wave vector, denoted in this 
paper by k^,, the single-particle spectral function involves 
a singular contribution at e = where e^, lies some 50 
to 70 meV below the Fermi energy e^, 

(ii) whether a uniform GS is sufficient, or inhomo- 

geneity, such as 'stripes' (whether static or dynamic) 
[63-65] are essential for the single-particle spectral func- 
tion Aa{k;s) to possess the properties indicated above 
and 

(iii) whether the above-mentioned properties expected 
of Acr{k; e) can arise through collective excitations inter- 
nal to the electronic system (charge and spin excitations) 
[66] or coupling of electrons to external bosonic modes 
(phonons) are essential, or a combination of both. 

As for the role that phonons can play in bringing about 

the above-mentioned properties concerning Aa{k:e), we 
believe that the experimental data considered by us in 
this paper [3-8] rule out a primary role; however, a sec- 
ondary (even though possibly vital) role is not excluded. 
As an example of the interesting possibilities that can 
become available through the mediation of the electron- 
phonon interaction, we have briefly discussed the Frohlich 
'shell' state. We have, however, shown that this state 
cannot be the appropriate state to which the experimen- 
tal observations as reported in [3-8] can be ascribed. 
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APPENDIX A: AN APPROXIMATE e^.^ 

In this appendix we deduce an approximate expression 
for e^.g. by assuming n^ik) e {0, 1}, subject to the condi- 
tion ria-{k) = No-, Vct, where N^- stands for the num- 
ber of particles with spin index a in the GS. The n^rik) 
pertaining to the (iVo- + -^CT)-pai"ticle GS of Hq in Eq. (37) 
is only a specific example of the type of ncr(fe) considered 
in this appendix. Many aspects in this appendix coincide 
with those in § A. 2 of appendix A in [19]; however, the 
approach adopted here is more direct and the details are 
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more closely related to the subject matters of the present 
paper. 

We consider the following form for the many-body in- 
teracting Hamiltonian H: 

a,cr' 

(Al) 

where ho{r) is the sum of the single-particle kinetic- 
energy operator and one-body external potential (e.g., 
ionic potential) and v{r — r') is the two-particle inter- 
action potential which wc leave unspecified at this stage. 
In [19] for the single-particle spectral function Aa-{r, r'\ e) 
pertaining to the GS of the Hamiltonian in Eq. (Al) we 
have obtained the following exact expression 

^ dee^^(r,r';£) = ha{r)V^^\r' u,ra) 

+ j d.'^r" v{r-r")^T^^\r'a,r"a';ra,r"a'), (A2) 

where r^'"^ m = 1, 2, is the m-particle ground-state den- 
sity matrix, for which we have (for details and various 
properties of this function we refer the reader to appendix 
B in [19]) 

r^™^(ri£7i, . . . , TmCTm; r[a[, r'^a'J 

(A3) 

the function T^^\r'a;ra) is more commonly denoted by 
Q<7{r',r). 

Now we specialize to cases in which one has the follow- 
ing spectral representation for the single-particle density 
matrix: 



(A4) 



where, for uniform GSs defined on the continuum, one 
has 



tpkir) 



1 



(A5) 



and, for uniform GSs defined on a regular lattice consist- 
ing of A^L sites (subject to periodic boundary condition) 



(A6) 



where {Rj || j = 1, . . . , N-s^} defines the above-mentioned 
lattice and (f>j{r) is a normalized 'atomic' orbital local- 
ized around r = 0. Since tpk{r) in Eq. (A6) is a Bloch 
function, it can be written as follows: 



(A7) 



where Uk (r) stands for a normalized periodic function of 
r satisfying Mfc(r) ~ Uk{r + Rj), Vj. In what follows we 
explicitly deal with the Bloch function in Eq. (A7) and 
deduce the results corresponding to the Bloch function 
in Eq. (A5) through identifying Uk{r) with unity. 
With 

ho{r)il>k{r) = Ekipkir), (A8) 
it follows that ^6 

ho{r) Bair', r) = ^ £fe n<,(fe) Mr) rk{r')- (A9) 

k 

For uniform GSs we have 

A,{r, r'; e) = ^ A,{k; e) MrWkir'), (AlO) 

k 

SO that through the orthogonality relation 

j d''rrk{r)M{r)=Sk,k', (All) 
from Eq. (AlO) we obtain 

A,{k;e) = j d^rd^r' rk{r)A„{r,r';e)Mr')- (A12) 

Thus from Eqs. (A2), (A9), (All) and (A12) we deduce 
that 

1 

-J deeA^{k;e)=ekna{k)+gl3<.^, (A13) 
where with reference to Eqs. (39) and (40) we have 
^ / d'rd^r' rkir)Mr') 
X y^dV w(r-r") X^r(2)(r'<T,r'V;ra,r'V), (A14) 



The RHS of Eq. (A2) is simply D^(r,r') defined in ap- 
pendix E of [19] (see in particular Eqs. (E5), (E12) and 
(E13) therein) of which the second contribution is equal 
to -B^{r,r') defined in appendix B of [19] (see Eq. (B29) 
therein). 



For time-reversal-symmetric GSs, for which £_f, no-(— fe) = 
£fc r\c{k) holds, one readily verifies that the RHS of Eq. (A9) 
is invariant under the exchange of r and r' and consequently 
so is also the LHS. 
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in which as in Eq. (37) wc have introduced the dimen- 
sionless two-body interaction potential w = v. 

From the representation in Eq. (A4) and the orthogo- 
naUty relation in Eq. (All) we have 

/ dV ga{r',r")gAr",r) = ^ {n^ik f ^rWkir'), 

k 

(A15) 

so that, in cases where ncr(fe) takes the values and 1 over 
the entire range of fe, it follows that ea{r',r) is idempo- 
tent (in the representation-free notation, idempotency is 
defined through g^g^ — g„). In such cases (and only 
in such case) the following result is exact: 



r^^^ (Xi,X2:,X-y,x'2) 



r(i)(a;i;a;'i) rW(a;i;4) 
r(i)(a;2;a;i) rW(a;2;4) 



(A16) 



in which Xi = Tiffi. By the conservation of spin we have 
[19] 

r(i)(r<7,rV) = 5,,,, r(i)(ra,r'<7) = 5,,,, ga{r,r'). 

(A17) 

Under the assumption that no-(fe) G {0, 1}, we thus have 
j dV w{r-r") ^r(2)(rV,r'V';ra,rV) 



-Vn{r; [n\)ga{r',r) 



w{r-r")Q,{r',r")Q,{r",r), 



(A18) 



where UH(r; [n]) is the Hartree potential corresponding 
to the total number density n(r) = '}2„n„{r), where 
no-(r) = Qaif-tr), that is 



^h(^; N) = / d'^r' v{r - r')n{r'). 



(A19) 



For two-body potentials for which J d'^r' v{r — r') = 
^v{q = 0) is unbounded, it is necessary to effect a reg- 
ularization procedure which takes account of the attrac- 
tive interaction of the particles with some uniformly dis- 
tributed background charge through which Vii{r\ [n]) is 



The fact that the single-particle density matrix pertaining 
to a non- interacting GS is idempotent is associated with the 
very specific condition where no-(fe) coincides with the char- 
acteristic function of the non-interacting Fermi sea (equal to 
unity inside and equal to zero outside the mentioned Fermi 
sea). The expressions that are presented in Eqs. (A20) and 
(A21) below are therefore of wider applicability than solely 
to the r'^' pertaining to non-interacting GSs. 



replaced by i'h(^'; [n — no]) where nQ:=N/n. For uniform 
GSs defined on the continuum, this procedure entirely 
removes Wh(^; [n]) from the formalism. For uniform GSs 
defined on a lattice, v^{r; [n]) is equal to a constant, in- 
dependent of r. Since in this appendix we deal with 
uniform GSs, in what follows we denote the Hartree po- 
tential by v-ii which represents a constant which may or 
may not be vanishing. 

Prom the above observations, some straightforward al- 
gebra, and Eqs. (A14) and (A18) we obtain 



/3^. = {^«H- y'^«)(||fc-fc'||)n.(fe')}n. 



Ef(fc)n4fe), 



(fc) 
(A20) 



where T,^^{k) stands for the Hartree- Fock part of the 

self-energy. For the uniform GSs of the single-band Hub- 
bard Hamiltonian, S^''(fc) = h~^Una (independent of k) 
where ria—Ng/Ni^; from this and Eq. (A20) it follows 
that /3<^ = nar)cr{k) and, through Eq. (40), £,k-M = rig. 
In arriving at Eq. (A20) we have considered the thermo- 
dynamic limit and replaced ^~^Y^ki by {2'k)~'^ j d'^fc'. 
This integral covers the entire available wave- vector space 
which for systems defined on a lattice consists of the IBZ 
corresponding to the underlying Bravais lattice; for these 
systems, the vector k — k' on the RHS of Eq. (A20) is 
to be understood as representing k — k' + Kq where the 
reciprocal lattice vector Kq (corresponding to an Umk- 
lapp process) is a function of k and k' and ensures that 
k~k' + KoG IBZ. 

From Eq. (A20) it is observed that within the frame- 
work of the approximation where na-(k) G {0,1}, f3'^.^ 
is proportional to ncr{k) so that for ^k;a as defined in 
Eq. (40) we have 



(A21) 



Thus, following Eq. (39), according to the present ap- 
proximation we have 



sk + h^fik). 



(A22) 



The unsound (see, however, next paragraph) nature of 
this expression in its functional form becomes evident by 
realizing the fact that by the assumed stability of the GS 
of the system it is required that e^.^. < /i, Vfe, while the 
expression on the RHS of Eq. (A22) in general does not 
satisfy this requirement. For instance, by considering 
the uniform electron-gas system for which one has the 
unbounded energy dispersion ek = ?i^||fc|p/[2me], one 
immediately observes that for sufficiently large ||fc|| the 
RHS of Eq. (A22) can be made to exceed any constant 
value. This aspect is made the more explicit by con- 
sidering the energy dispersion e^.^ in Eq. (48) (sec also 
Eq. (23)). Under the conditions for which Eq. (A20) ap- 
plies, Eq. (48) yields 



25 



(A23) 



References 



which identically coincides with the (approximate) result 
in Eq. (A22). This result, similar to that in Eq. (A22), 
clearly exposes the shortcoming of the approximation ac- 
cording to which r\^{k) G {0, 1}; in [9,10] is has been 
shown that the exact s^.^ satisfies s^.^ > ji, Vfe. 

The peculiarities of the expressions in Eqs. (A22) and 
(A23) arc understood by considering the fictitious single- 
particle spectral function Aa{k; e) introduced in Eq. (32) 
and to which the energies e^.^ and e^.^ correspond. It is 
observed that Aa{k:e) yields the same result by substi- 
tuting Aa{k;e) for the exact A„{k;e) in Eq. (39), that 



IS 



/^^de £^^(fc;e) 



Vfe. 



(A24) 



The same applies to e^.^ in Eq. (23). As is evident, the 
validity of Eq. (A24) and of its counterpart concerning 
^fe o- crucially depends on the properties e^.^ < fi and 
e^.^ > ;U, Vfe. In order to gain insight into the conse- 
quences of the approximate framework considered in this 
appendix, leading to identical energy dispersions for e^.^ 
and £fc.CT, for all fc, it is appropriate to substitute the 
expressions in Eqs. (A22) and (A23) into Eq. (32) from 
which it follows that 



where 



Aa{k;e)=n5{e-ef.„), 



(A25) 



(A26) 



One observes that, within the approximate framework in 
which n^ik) € {0,1}, the single-particle energy disper- 
sion in fact consists of a single branch, namely £^.^, and 
not of two identical branches. Thus the identity of e^.^ 
and e^,^ in our approximate framework amounts to no 
fundamental defect, rather to a manifestation of the re- 
stricted space of the single-particle excitations available 
to states (or GSs) for which the associated ^^{k) satis- 
fies ^^{k) G {0, 1}, Vfe. With reference to specifically our 
considerations in § 4, resulting in the conclusion that for 

cases in which v{r — r') is the long-range Coulomb poten- 

< 

tial (or one as long-ranged as this), Vfe£^.^ is logarith- 
mically divergent for fe approaching regions (e.g. S^-a) in 
which ncr(fc) is discontinuous, it is observed that this as- 
pect is appropriately preserved by ^k^^-u in which e^f^ 

is purported to approximate s^.^. □ 
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FIG. 1. Schematic representation of the conditions under 
which the quasi-particle equation + h'Ecr{k;e) = e has ap- 
parently two solutions at fc = fc*, where fe* is located strictly 
inside the Fermi sea FSa', the two solutions are denoted by e^- .^ 
and £^+.^, where e^.-.^ < M £^+.^ > M- The singularity 
of n„{k) at fe = fe, together with the normalization condition 
/_oo Aa{k; e) = 1 imply that e^- and £^+.^ must be 
solutions of the quasi-particle equation at fe = fej and fe — kt 
respectively (see § 3.2). Curve (a) depicts e^,^ +h'Re[T,cr{kf ; e)], 
curve (b) gives hlm[T.c{kf;e)], and curve (c) shows the func- 
tion f{e):=e. Although the functions are depicted schemati- 
cally, care has been taken to ensure that film\Tia{k^;e)] sat- 
isfy the necessary requirements Im[Eo-(fef ; e)] > for e < /x, 
Im[E<^(feJ;e)] < for e > ^, and Im[S„(fef ; /i)] = 
(the stability of the GS requires the more stringent condition 
Im[ECT(fe; /i)] = 0, Vfe). Since nCT(fe) is assumed to be discon- 
tinuous at fe = fe*, Tia{k^;e) is continuously differentiable with 
respect to e in the neighbourhood of e = ^k^-a- that, 
since following the considerations in § 3.2 for the case dealt 
with here, the spectral weights Z^-,^ and ^^+.^ correspond- 
ing to the singular contributions SCT(e) and respectively 
are equal (we denote the common value by Z/^^.^), it follows 
that dECT(fe7 ; e) /de at e = ^fc- equal to dEo-(fe^; e)/de at 
e = Eig+.^- Further, since in general < Z)^^.^ < 1, it follows 
that the latter derivatives are negative or at most zero (the latter 
corresponding to Z)^^.^ = 1). 
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FIG. 2. Schematic representation of Ej^.^, the asymptotic so- 
lution of the quasi-particle equation (Eq. (3)), and as de- 
fined in Eq. (39), for fe inside the Fermi sea and close to the 
Fermi wave vector fepio-, with fepio- — fe pointing in the direc- 
tion of the outward normal n{k^'a) to the Fermi surface Sf-ct at 
fe = kp-cr- It is assumed that nCT(fe) is continuous at fe = k^^cr 
(i.e. ^fcF.„ = 0) while it undergoes a finite discontinuity at 
k — ki, (i.e. Z]^^.„ > 0). Following Eqs. (4) and (27), e^.^ and 
e^.^ coincide at fe = kY-<y and k = ki, respectively and, owing 
to the smallness of ||fe, — feF;CT||/||feF;(T|| (in practice, of the order 
of 5%), the possible difference between e^,^ and cannot be 
experimentally discernible for fe in the interval between fe* and 
fep^CT. We have ^ = tan^^ (7 rj,(feF;c7) ■ VfcE^.^I ), in which 

7 = ||feF:CT||/£F, and (j> ~ —Zf^^.^9 /na{k^); counterclockwise is 
the direction in which we count angles as positive. By viewing 
£^.^ as the single-particle energy dispersion as measured through 
ARPES and excluding the results corresponding to the under- 
doped cuprates, our schematic representation (which is based on 
Eq. (63); see also Eqs. (74) and (75)) is in excellent qualitative 
agreement with the experimental observations by Zhou et al. 
[3]. For fermions interacting through the Coulomb interaction 
potential, Ti(feF;a) • ^k^k-^ scquires a logarithmic divergence at 
k — ki, (see Fig. 3). Taking this aspect into account, as well as 
a finite discontinuity in £^.^ at fe = fe,, which is feasible (§ 3.4), 
the corresponding schematic representation would similarly be in 
excellent (qualitative) agreement with the experimental observa- 
tions by Zhou ct al. [3] concerning cuprates in the underdoped 
regime, as exemplified by (La2_a;Sr3;)Cu04 for x — 0.03. 
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FIG. 3. The energy dispersion e^.^ = Sk + g(,k:^ (Eq. (39)) 
for n:=N/Ni, equal to 1 (curve (a), corresponding to the ex- 
treme limit of hole underdoping) and 0.7 (curve (b), corre- 
sponding to an overdoped cuprate compound) along a diago- 
nal direction of a square IBZ as calculated within the frame- 
work in which g^fe;CT is replaced by ?iE™(fe) (for the details 
concerning this choice see the § 4.1 and appendix A; see also 
Fig. 2) evaluated in terms of the long-range Coulomb potential 
and a GS momentum distribution function ncr{k) of the form 
ncr(xkY-cr) = 1 — hQQ{x — xq) — Zy.^,^Q(x — 1) where < a;o < 1 
and ho = (1 - ■^fcp.„)(a;z ~ '^)l'yx\ - x\); here fez = izk^;^ 
stands for the vector from the centre of the IBZ to the zone 
boundary in the direction of kp-.a G 5f;ct. Full account has been 
taken of the Umklapp of wave vectors that in the process of 
wave-vector integration move outside the IBZ. The calculations 
are performed for xq = 0.95 and Z^^.^ ~ 0, VfcFio-, so that by 
normalization ho = 0.968 for n = 1 and ho = 0.976 for n = 0.7. 
For the energy dispersions (curve (c) for n = 1 and curve 
(d) for n — 0.7) we have chosen the tight-binding expression 
£^ = —2t[cos{aokx) + cos{aoky)] where kx and ky stand for the 
Cartesian coordinates of fe and ao — 3.895 x 10~^" m for the 
lattice constant. For the Fourier transform of the Coulomb po- 
tential we have wdlqH) = giu(||(7||) where it)(||(jr||) = 27rao/||<j'| 
and g — / (Aiveo^rao) in which stands for the square of 
the electron charge, eo = 8.854... X 10"^^ F m~i for the vac- 
uum permittivity and the relative dielectric constant of the 
background. With eo:=?i^/(m*ao), the calculations presented 
here correspond to g = 3eo. Here m* denotes the effective 
mass associated with (to be distinguished from that associ- 
ated with e^.jj.) at fc = fepiCT along a diagonal direction of the 
IBZ, for which we have chosen m* — ArUe (for both densities 
n — 1 and n = 0.7), where me stands for the bare electron 
mass (for the above values of ao and m*, eo = 0.1255 eV; 
the choice g = 3eo thus amounts to = 9.814). Consequently, 
t = eoKF/[2\/2 sin(KF/\/2)], where KF:=ao||feF;CT||. For n = 1 we 
have Kf = 2.221 {hip/[V2n] = 0.5) and for n = 0.7, kf = 1.965 
[k^/IV^tt] = 0.442), corresponding to t ^ 0.785eo for n = 1 
and t — 0.706eo for n = 0.7. A density-dependent t is chosen 
to avoid a strongly density-dependent m*. 
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FIG. 4^ The analytic structure of the single-particle Green 
function Ga{k;z) on the physical Riemann sheet of the com- 
plex z plane for cases where (A) fc — k^.^., with feF;o- G 5f;ct 
and (B) fc G FSa and at a finite non-vanishing distance from 
5f;ct. It is assumed that in the latter case for the fc at is- 
sue the quasi-particle equation + hT,„{k;e) = e is satisfied 
for (fc,e) = (fc'^,efc=F;<,) (see Fig. 1) in which e^T-^ = ej.^ 
(Eq. (13)); for metallic states and fcF;o- G 5f;ct, this equation is 
satisfied for (fe,e) = {k^.^^,eY)- The solid lines on the real axes 
whose end points coincide with /i^.^:=i?jv^,jv^;o — -Ejv„-i,jv^;o 
and j-L^.^:=EN„+i,Nff;0 — EN„.Ng;0, represent the branch cuts of 
Ga(k;z); here Em^^Mh-.o, with M„ = N^, + stands 

for the total energy of the (A/o- -I- AfCT)-particle GS of the interact- 
ing Hamiltonian H (for a a —l, and vice versa). For metallic 
GSs, /i^.^ = Ef and ^t^.^ differs only infinitesimally from M^.^; 
the chemical potential /i satisfies A^iv-o- < A* < f^tr-a-- evaluat- 
ing the integrals of z"^G„{k; z), m = 0, 1, over the contour C 
we employ the decomposition dz = Jq\^q dz + dz which 
amounts to the decomposition, according to Eq. (8), of Acr{k; e) 
into its regular and singular contributions respectively. 
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